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SUMMARY 


Typical pollutant gas concentrations at the stack exits of 
stationary sources can be estimated to be about 500 ppm under 
the present emission standards. Raman lidar has a number of 
advantages which makes it a valuable tool for remote measure- 
ments of these stack emissions. Tests of the Langley Research 
Center Raman lidar at a calibration tank indicate that night 
measurements' of SO 2 concentrations and stack opacity are possible 
Accuracies of 10 percent are shown to be achievable from a dis- 
tance of 300 m within 30 min integration times for 500 ppm SO 2 
at the stack exits. All possible interferences were examined 
quantitatively (except for the fluorescence of aerosols in actual 
stack emissions) and found to have negligible effect on the 
measurements. An early test at an instrumented stack is 
strongly recommended. 
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REMOTE MEASUREMENTS OF POLLUTANTS FROM STATIONARY 
SOURCES USING RAMAN LIDAR 

By 

S.K. Poultney^ , M.L. Brumfield^, 
and J.S. Siviter^ 

1 . INTRODUCTION 

Two types of stationary sources contribute to regional 
pollution throughout the country if adequate controls are not 
used. These are fossil fuel power plants and municipal inciner- 
ators. Smelters, pulp mills, and petroleum refineries can also 
contribute in particular locations. The emissions which lead to 
pollution consist of both gases and particulates. The primary 
gases are S02» NO, and HCl. Subsequent dispersion and reaction 
of these components and others can lead to sulfurous smogs and 
contribute to photochemical smogs with their attendant effects 
on man and property. Regional air quality standards have led to 
emission standards for existing stationary sources. Estimates 
of the emissions from the various sources can be made by knowing 
general information about the activity at the source, the type 
of fuel, and the control devices in use. Many of the operators 
of stationary sources will monitor the emissions to stay within 
emission standards and to check the efficacy of their control 
devices. However, the general community may want to verify the 
compliance of each stationary source with the emission standards 
by independent measures from a remote location. 


^ Department of Physics and Geophysical Sciences, Old Dominion 
University, Norfolk, Virginia 23508. 

NASA Langley Research Center, Hampton, Virginia 23665. 
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Raman lidar has a number of advantages which makes it a 
valuable tool for remote measurements of these stack emissions. 

The typical pollutant gas concentrations at the stack exit 
(e,g./ 500 ppm in the absence of controls) can be monitored by 
means of Raman- scattered light from the illuminating laser beam 
in spite of the rather low cross sections for the scattering 
process. The measurements are possible at the typical stack 
height of several hundred meters from a remote location and are 
specific to the gas in question. The illuminating laser need not 
be tuned nor narrowed as in other schemes. Supplementary measure- 
me^nts of stack plume cross section and stack exit velocity must be 
made in order to derive the concentration and mass emission rate. 
The plume diameter can be obtained using a separate channel of 
the Raman lidar. Lidar calibration to obtain absolute values can 
be achieved in several ways . One can operate two or more simul- 
taneous return channels with one being ambient nitrogen outside 
the plume. The relative sensitivity of the channels can be 
periodically checked by viewing a standard lamp. Altern^ately , 
one could operate a single channel viewing the pollutant in the 
plxime and then oxygen outside the plume. This latter method is 
most convenient for SO2 measurements. 

Two types of interference must be guarded against : leakage 

of the wrong wavelength return or acceptance of aerosol broad-band 
fluorescence into the spectral bandwidth of a receiver channel. 

The leakage can be controlled by proper design of the spectral 
selection elements in the lidar receiver. Fluorescence may not 
be a problem for stacks with particulate control devices or for 
longer wavelength lasers. It can be identified by looking for 
scattered light at wavelengths not characteristic of Raman scatter 
and for signals decaying with a characteristic fluorescence time 
constant. Both of these possible interferences as well as the 
actual performance of the Raman stack monitor lidar can be checked 
at a calibration tank facility before measurements at a stack. 

This report describes a particular Raman lidar and its performance 
at the calibration tank facility. This Raman lidar is capable of 


night measurements of SO 2 and HCl under the conditions mentioned 
above to an accuracy of 10 percent within 30 min from a distance 
of 300 m. Its performance for SO 2 has been independently measured 
in a calibration facility for the first time. This report does 
not address directly the issue of what laser should be used for 
the optimum Raman lidar for either performance or eye safety. 

The particular lidar described here can be used for day measure- 
ments only if a different data acquisition scheme were to be used. 

The measurement of stack particulates most accessible to a 
Raman lidar is that of stack opacity. The Raman return from 
nitrogen (or oxygen) in the range cells before the stack and the 
range cells after the stack can yield the extinction or opacity 
of the stack plume. The use of a Raman lidar for this measurement 
may eliminate certain problems associated with other lidar means . 
Concentration and mass emission rates are more difficult to obtain 
by lidar means. In addition to the supplementary measurements 
needed for pollutants, one must know or measure the optical prop- 
erties of the particles involved. Lidar can easily provide only 
two parameters related to particle concentration and properties; 
extinction and backscatter. The measurement of particulate mass 
emission rate has been addressed by several authors, but will not 
be discussed here. 

2. STATIONARY SOURCE EMISSIONS AND EMISSION STANDARDS 

2A. Estimates of Emissions from Power Plants 
and Municipal Incinerators 

Emissions from the stationary source fossil fuel power plants 
and municipal incinerators will be discussed in this section. The 
amount of these emissions can be estimated by knowing general infor- 
mation about the activity, fuel, and control devices of the par- 
ticular source (ref. 1). Uncertainty about any of these components 
or their use requires remote measurements to check independently 
compliance with emission standards. Fossil fuel power plants use 
gas, oil, coal, or a combination of these three fuels. A municipal 
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incinerator might also be used as a power plant. Its fuel would 
be rubbish. Each fuel has an emission factor which indicates 
how much pollutant is emitted for a given amount of fuel burned. 

This factor in the case of SO 2 depends on the sulfur content of 
the particular lot of fuel. Table 1 lists the emission factors 
and caloric values for the several fuels except gas which is very 
clean. The percent sulfur, fly ash, and chlorine in the respective 
entries is represented by the letter s. For example, 1 percent 
sulfur coal would yield 2.7 kg of SO 2 for the conditions stated 
in table 1, The amount of N0„ emitted depends on the flame 
temperature at which the fuel burns. Table 1 also lists the 
actual amount of pollutant emitted per 4.186 GJ heat input (also 
written GJh) . The emitted amount in grams will be seen to be 
closely related to the present emission standards. The above 
heat input each second represents a large modern 1.7 GWe power 
plant, assuming an efficiency of 40 percent. 

One can estimate the pollutant mass emission rates and con- 
centrations near the top of the stack if one knows the rate at 
which the power plant in question uses fuel. Table 2 lists the 
rates and concentrations for the modern 1.7 GWe power plant. 
Calculation of the mass emission rates is straightforward. Cal- 
culation of the concentration requires additional assumptions. 

The concentration in the plume at the stack exit is given by 

C = Q^/vA (1) 

where C is the concentration in ppm, the volume emission 

rate in m^/sec, v the exit velocity in m/sec, and A the plume 
area in m^ right above the stack exit. It is assumed here that 
V equals 100 km/hr and A corresponds to a stack diameter of 
10 m. The volume emission rates for gases are obtained from the 
mass emission rates using 

= 22.4 x 10“3 (Q^/Mp) (T/273) (2) 
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Table 1. Stationary Source Emissions. Values taken or calculated from Williamson 
(ref, 1). The sulfur, chlorine, and fly ash content of actual fuel used 
is specified by s in percent. The amount of pollutant emitted is cal- 
culated for a heat input of 4.186 GJ, which is typical of modern large 
power plants with 1.7 GWe capacities. 



Petroleum 

Coal 

Rubbish 

Emission Factors 

SO 2 

19s g/liter 

19s kg/tonne 

1 kg/tonne 

NO 2 

12 g/liter 

10 kg/ tonne 

1 kg/tonne 

Particulates 

1.2 g/liter 

8s kg/ tonne 

15 kg/tonne 

HCl 

— 

10s kg/ tonne 

10s kg/tonne 

Caloric Value. 

4.2 X 10^ J/liter 

30 X 10® J/tonne 

11 X 10® J/tonne 

Amount of Pollutant 
Emitted 

SO 2 

2.0s kg 

2.7s kg 

0.38 kg 

NO 2 

0.13 kg 

1 . 4 kg 

0.38 kg 

Particulates 

0*12 kg 

8.2s kg 

5.7 kg 

HGl 


1.3s kg 

3.8s kg 



■ f ■■ “ 

a 




Table 

2. Fossil Fuel Power Plant Mass Emission Rates and Plume Pollutant 
Concentrations at Stack Exit. Values calculated from Williamso: 
(ref. 1). Plant is 1.7 GWe capacity operating at 40 perceht 
efficiency with a 10 m diameter plume at exit and a 100 km/hr 
exit velocity. Stacks are assumed to have a 99*7 percent par- 
ticle precipitator. 


SO 2 

NO 

X 

Particulates 

HCl 

Fuel 





Petroleum 

400 ppm 

38 ppm 

— 

— 

1% S 

■ (2.0 kg/sec) 

(0.13 kg/sec) 

(0.36 g/sec) 

— 

Coal 





1% S 

540 ppm 

380 ppm 

— 

230 ppm 

1% fly ash 

(2.7 kg/sec) 

(1.4 kg/sec) 

(49 g/sec) 

(0.68 kg/sec) 

0.5% Cl 





Rubbish 





3% plastic 

7 5 ppm 

110 ppm 

— 

680 ppm 

0.5% Cl 

(0.38 kg/sec) 

(0.38 kg/sec) 

(17 g/sec) 

(1.9 kg/sec) 








are 


where Q indicates the respective rates, the units of Q„ 

M 

kg/sec, Mp is the molecular weight of the pollutant in kg, 
and T is the exit plume temperature in ®K. It is assumed here 
that the plume temperature is 413® K. The fuel is assumed to be 
1 percent sulfur and, in the case of coal, 1 percent fly ash. 
Similar calculations are done for HCl from plants burning 0.5 
percent Cl coal and rubbish with 3 percent plastic component. 

Note that the HCl emission from a municipal incinerator may be 
an order of magnitude higher if plastics with high polyvinyl 
chloride components are burned. 

The following observations can be made based on table 2. 

The SO 2 concentrations at the top of large, modern power plants 
(e.g., 1.7 GWe) without controls are about 500 ppm for 1 percent 
sulfur fuel which is just in the range of remote Raman lidar. 

The difficulty of disposing of acid liquid wastes means that 
these SO 2 concentrations will not be reduced by controls except 
if the sulfur is removed from the fuel before use or if it can 
be commercially recovered from the wastes right at the plant. 
Relaxation of present emission standards for coal plants due to 
the energy crisis could increase SO 2 concentrations dramatically 
since coal can have sulfur contents as high as 7 percent. Such 
relaxation would, however, require either modification of stack 
design or use of control devices because the increased emissions 
may impinge on the ambient air quality standards as discussed in 
section 2C. The NO^ concentration is relatively high for a coal 
stack, but it turns out that the Raman cross section for NO (the 
form of NO at the stack exit) is an order of magnitude below 
that of SO 2 . Raman lidar measurements of NO are therefore not 
promising. Table 2 indicates that HCl up to the 250 ppm level 
might be present in coal stack emissions for 0.5 percent Cl coal. 
The higher HCl concentrations (e.g., 700 ppm) in the plume from 
a municipal incinerator power plant of the dimensions stated 
above should be even more accessible to the Raman lidar. Particu 
late mass emission rates are also given under the assumption of 
99.7 percent efficient precipitators. 



The above estimates for stack plume concentrations of pollu- 
tants have been made for large modern power plants operating at 
full capacity. One might ask what concentrations are expected 
for the older 200 MWe power plants and what concentrations are 
expected when either size plant reduces its power output. By 
equation (1) , the expected concentration depends on the mass 
emission rate, the exit velocity, and the stack area. The mass 
emission rate has an obvious relation to power plant capacity and 
output. Changes in operating power levels at a plant will change 
the concentrations of pollutants in the plume directly with the 
change in emission rate. if the exit velocity is kept the same. 

At many plants, the exit velocity is decreased as the output 
level is decreased so that the pollutant concentration in the 
plume may stay the same. A measurement of stack exit velocity 
is essential if Raman lidar concentration measurements are to 
be related to the required mass emission rates. The 200 MWe 
power plants have one-tenth the mass emission rates of the large 
plants by definition, but their stack sizes are scaled down from 
the sizes of the stacks at the large plants for reasons discussed 
in section 2C. Typical stack diameters of the small plants are 
one-third those of the large plants. Equation (1) thus indicates 
that the pollutant concentrations are the same for large and small 
plants which use the same fuels. A measurement of stack plume 
cross section is essential to obtain the required mass emission 
rates from the Raman lidar concentration measurements (in addition 
to interpreting the lidar measurements in the first place) . 

2B. Emission Standards 

The Federal standard of performance for fossil-fuel-fired 
steam generators with respect to SO 2 emission is quoted in grams 
of SO 2 in the emitted gases per million calorie heat input (Mcalh) 
(ref. 2). For liquid fossil fuel, the number is 1.4 g/Mcalh, 
and for solid fossil fuel, the number is 2.2 g/Mcalh. Table 1 
lists the pollutant amount emitted per 1 Gcal heat input or 
4,2 GJh for various fuels. To compare these estimates with the 



Federal emission standards, one merely divides the amounts in 
table 1 by 10 For example, table 1 indicates that 1 percent 
sulfur coal emits slightly more SO 2 than allowed by the emission 
standard* The Federal emission standards are apparently a restric- 
tion on the sulfur content of the fuel used. The standards also 
do not limit the absolute amount of pollutants which a plant may 
emit since a plant with a larger power capacity requires a larger 
heat input per second. The emission standards are, of course, 
derived from ambient air quality considerations and do not contain 
the final word on stack emissions. They are a guide to the achieve- 
ment of ambient air quality standards and as such are contested 
for various reasons by stationary source operators. It will be 
useful to briefly consider the issues involved to see what role 
Raman lidar may play in resolving the emission standards contro- 
versy. 

2C. Dispersal of Stationary Source Emissions and 
Ambient Air Quality Standards 

The ambient air quality standards of 1971 (ref. 1) are summar- 
ized in table 3 for selected pollutants. Little controversy 
surrounds these standards . Note that the standards are concen- 
trations in ppm and are much lower than the stack concentrations 
(e.g., 0.1 ppm for SO 2 ). Plume dispersion processes are thus 
very important. Conventional plume dispersion theory predjicts 
that the maximum ground-level pollution concentration depends 
directly on the pollutant emission rate (not stack concentration) 
and inversely on the square of the stack height (ref. 1) . The 
increase in mass emission rate allowed the larger capacity power 
plants can thus be offset by proper stack design. For example, 
the 200 MWe plants have stacks about 100 m high, whereas the new 
1.7 GWe plants have stacks about 300 m high. The maximum ground 
level pollutant concentration of each is probably about the 
same. The actual value of the ground-level concentration may 
vary widely. The dispersal process depends on the geographical 
location of the plant, the meteorological parameters of the 
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Table 3. 


Selected Air Quality Standards (Federal) 1971, From Williamson (ref, 1) 


Pollutant 

SO2 

NO2 


Averaging Time 
Annual Average 
2 4 -hour 

Annual Average 


Primary Standard 
0.03 ppm 
0.14 ppm 
0.05 ppm 


Suspended 

Particulate 

Matter 


Annual Geometric 
Mean 

2 4-hour 


75 microgram/m^ 
260 micrograra/m^ 


Objectives of Standards 

To prevent possible increase 
in chronic respiratory disease 
and damage to vegetation. 

To prevent possible risk to 
public health and atmos- 
pheric discoloration. 

To improve visibility and 
prevent acute illness when 
present with about 0.05 ppm 
SO2 . 



atmosphere (including turbulence ) , the design of the plant and 
the stack, and chemical reactions (in the presence of particu- 
lates) that may be taking place. Even in the best of conditions, 
one does not expect predictions accurate to better than a factor 
of two over a certain averaging period. Higher concentrations 
over shorter times at the ground could be dangerous and could 
occur. Williamson (ref. 1) concludes his tutorial review of 
plume dispersion as follows; "A useful quantitative theory per- 
mitting the exact calculation of the downwind pollution from a 
point source, which incorporates measurable features of atmos- 
pheric turbulence in the intervening distance, has not yet been 
achieved.” 

Federal emission standards and designs for new power plants 
were undoubtedly determined by application of the present dis- 
persion theory and by extrapolation from empirical measurements 
in the vicinity of existing plants. Either means raises questions 
which lead to the present controversy about emissions standards. 
First, how much credence can one place in emission standards 
derived from ambient air quality standards via an imperfect theory 
Part of the problem in constructing a better theory is that the 
rapid temporal variations of the ambient concentration of pollu- 
tants cannot even be precisely measured. Neither can the meteoro- 
logical parameters be measured over the time scale and over the 
affected area. In view of these measurement difficulties, it 
is natural to ask whether or not Raman lidar could obtain the 
necessary accuracy, range, and time resolution. For 1 ppm SO 2 , 
the use of a Raman lidar such as the one .described below would 
require a decrease in range to about 10 m, thereby making the 
lidar no longer remote. The same conclusion is true for the 
several meteorological parameters accessible to Raman lidar. 

The empirical means for determining the emission standard 
immediately raises the question of why a standard other than the 
ambient air quality standard. Each point source would maintain 
a monitoring network and would operate so as to never broach the 
ambient standard. On good days, the source would burn high sulfur 


fuel, and at bad times it would switch to low sulfur fuels (ref. 3). 
The problem with this approach is that pollution is usually caused 
by a group of point sources. The present emission standard at 
least means that each individual in this group maintain its pro- 
portionate share of pollution. In fact, a variable emission 
standard based on atmospheric conditions might not be a bad 
compromise. In any case, one might like to have an instrument 
that can check for compliance the operation of each individual 
point source. The Raman lidar is such an instrument. 

A third possible use of Raman lidar for stack plume monitoring 
might come in the chemical reaction aspect of the dispersion process. 
For example, the SO 2 concentration might be monitored in the presence 
of varying amounts of particulates or other material. The lidar 
might, at the same time, monitor the particulate matter present. 

3. RAMAN LIDAR 

Raman lidar is not a new idea. Its use was first suggested 
by Cooney (ref. 4) in 1965 and implemented for the detection of 
atmospheric nitrogen by Leonard (ref. 5) in 1967. Its applica- 
tion to the measurement of pollutants was first proposed by 
Inaba and Kobayasi (ref. 6) and then realized by them (ref. 7) 
in 1970 for the cases of SO 2 and CO 2 in a nearby oil smoke plume. 

A brief summary of the use of Raman lidar for stack emission 
measurements is given below. First, however, the Raman scattering 
phenomenon will be described. Then follow a description of the 
LaRC Raman lidar used in the present tests and predictions of 
its performance for various gases. After the lidar summary, the 
issues of possible interferences and absolute calibration are 
outlined. 

3A. The Raman Scattering Phenomenon 

Raman scattering may be regarded as an inelastic collision 
between a photon of energy hv^ and a molecule. The molecule 
is left in an excited state characteristic of its structure and 
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a photon with energy hvj^ = hv^ ” AE is emitted. The energy 
loss AE of the photon corresponds to a unique rotational or 
a vibration-rotation energy change in the state of the molecule. 

The process is sketched in figure 1 for vibrational Raman scatter- 
ing. In the wavelentjch domain, this shift is to the long wave- 
length side of the exciting laser wavelength as indicated in 
figure 2. This paper concerns itself with the large shifts of 
the vibration-rotation excitation. Figure 3 indicates specific 
shifts for the constituents of an oil smoke plume for excitation 
by a ruby laser. Table 4 lists these shifts, the scattered wave- 
lengths for ruby excitation, and more accurate scattering cross 
sections for a number of gases of possible interest. The cross 
sections have been stated relative to that of nitrogen and the 
appropriate corrections have been made for the wavelength depend- 
ence of Raman scattering using equation (3) and the actual exci- 
tation wavelength of measurement. 

a ~ 1/X^ (3) 

These cross sections are quite small. For comparison, the Rayleigh 
scattering cross section is 1.98 x lO"^® cm^/sr (ref. 11). The 
Raman cross sections are the prime determinate of the design and 
performance of a Raman lidar. 

3B. Raman Lidar Design 

The design of Raman lidars is not particularly straightforward 
due in part to the weak scattering cross sections and in part 
to many possible sources of interference. Kildal and Byer (ref. 12) 
and the work of the Japanese group (refs. 13, 14) give a good 
insight into most of the intricacies of the design although each 
omits certain possible interferences. A Raman lidar consists of 
a laser to illuminate the target gas, a telescope receiver to 
collect and separate the scattered light, and a data acquisition 
system (DAQ) to record the data. Figure 4 shows the first two 
components schematically and figure 5 shows the DAQ. The choice 
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Table 4. Raman Scatter Properties of Selected^ Gases. Q branch cross sections are 
given relative to nitrogen for 6943 A excitation. The appropriate wave- 
length scaling was done from the wavelengths of actual measurement. 


Gas 


Raman Shift 
(cm~^) 

(6943) 

(A) 

(6943) 

Reference 

N2(Vj) 


2331 

8284 

1.0 

8 

02 


1556 

7784 

1.4 

8 

SO2 


1151 

7546 

6.6 

8 

NO 


1876 

7982 

0.51 

8 

CO 2 (2V2 

) 

1286 

7624 

1.10 

8 

CO 2 (Vi ) 


1388 

7683 

1.59 

8 

CO 


2143 

8158 

1.02 

8 

H 2 S 


2611 

8480 

6.4 

9 

CH^Cvi) 


2914 

8703 

7.5 

9 

HCl 


2886 

8683 

2.6 

10 

Note : 

Absolute N 2 differential scattering 

cross section is 

• 

• 


4.4 

(± 1.7) X 10“3i 

cm^/sr for 5145 

0 

A excitation (9) 



1.1 

X 10"^^ cm^/sr 

0 

for 6943 A excitation. 





of components is dictated to a large extent by the lidar equation 
(4) which predicts the signal return from a target gas based on 
lidar system parameters and on the concentration and cross-section 
of the gas. Elimination of interferences is addressed in section 
5E. 


The lidar equation is here expressed in terms appropriate 
for photon counting of the return signal. It is divided into 
four dimensionless factors. 



First, there is a scale factor which adjusts the conventional 
lidar equation (usually expressed in terms of instantaneous 
power) so that the scattered signal can be expressed in counts 
per range time bin and the exciting pulse in photons. The second 
factor is the number of photons in the exciting light pulse of 
E joules leaving the laser. The third factor contains the 
backscatter cross section, 02 (it) , the gas concentration, n 2 , 
the length of the viewed interval, L, and any extinction along 
the path due to aerosol, 3a(tot), or absorbing gas n a . 

This factor assumes that the viewed interval is much longer than 
the exciting pulse length. Corrections are necessary if the 
range interval is shorter, depending on whether the gas is 
localized in a plume or continuously distributed. The last 
factor contains the total efficiency of a receiver channel, 

K, including the photosurface quantum efficiency and it con- 
tains the lidar geometry with A the area of the receiver and 
R the range to the viewed interval. 

The lidar parameters used or assumed in the present design 
are listed in table 5. The total channel efficiency is given 
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Table 5 


LaRC Raman Lidar Design Parameters 


Range 

Range interval 

SO 2 concentration 

Laser wavelength 

Laser energy 

Laser beam divergence 

Laser firing rate 

Receiver diameter 

Receiver field of view 

Detection method 

SO 2 Raman channel total 
efficiency, K 

Atmospheric transmission 
(two way) 


300 m 

6 m or 40 nsec 
1000 ppm 
6943 A 

1.5 J or 6 X 10^® photons (nominal) 
1 mrad 
30 ppmin 

8 inches (or 0.03 m^ in area) 

3 mrad 

Photon counting 

0.025 cts/scattered photon 

0.8 


only for the most sensitive channel viewing SO 2 . The design is 
the result of a number of separate contributions which# however# 
had the same basic aims. The smallest telescope possible was 
chosen for compactness. This choice both required and allowed 
photon counting techniques to be used. Use of photon counting 
detection also allowed convenient use of digital elements in 
the DAQ (ref. 15) and narrow range bins'(e.g.# down to 20 nsec) 
(ref. 15), but limited lidar operation to night-time as will be 
discussed later. The DAQ configuration as shown in figure 5 
includes the Tomlinson Research sixteen-bin photon counting unit 
(ref. 15) and an auxiliary single time bin second channel which 
was used for a time. The spectral selection elements were chosen 
to be interference filters as shown in figure 4 for compactness, 
ease of operation# and higher throughput than many high-rejection- 
ratio monochromators. The selection of a ruby laser gave the 

O O 

option of operation at 6943 A or at 3472 A at reduced energy 

O 

per pulse. The choice here was 6943 A due to the simplicity of 
operation, the higher pulse energies available# and the hopes 
of avoiding various interference problems. 

The original version of the LaRC Raman lidar was described 
by Melfi et al. (ref. 16) and is basically the same as herein 
described except for certain changes and improvements. Northam 
and Brumfield (ref. 17) later reassembled the lidar at LaRG 
using a more compact 8-inch telescope because of the ND 1.0 
attenuator found necessary by Melfi with the original 24-inch 
receiver. The change in telescope required Northam and Brumfield 
to add a new detector package. Northam also initiated plans to 
evaluate the lidar in the calibration tank mode. Siviter assembled 
the calibration tank. The lead author then made contributions 
to the reliability of photon counting# the efficiency of detection, 
the certainty of specific pollutant detection, the flexibility of 
the data acquisition system# formulation of data analysis programs, 
and direction of the tests reported here. The photon counting 
and other detector contributions are outlined in Appendix 2 . 
Efficiency of detection was improved by addition of the high 


17 


efficiency RCA 31034 photomultipliers, the acquiring of new, 
hybrid filter elements, and by the addition of provisions to 
tilt the interference filters to their peak transmission. The 
filters were carefully specified to be centered slightly above 
their Raman wavelengths so that the tuning distance would be 
small. Certainty of detection was improved by means of this 
ability to tilt the filters and by means of the hybrid filter 
combinations. The flexibility of the DAQ was improved by a 
hardware change to allow a dual (but interleaved) channel mode 
(ref. 15) . Diagnostic and data analysis programs were formulated 
and appear in Appendix 1. The results to be presented in section 
5 were directed by the lead author. Several changes in the LaRC 
lidar would have improved its performance still further, but time 
was not available to do both the demanding evaluation of the 
lidar at the calibration tank and the changes. These changes 
include an analog DAQ for daytime operation, better beam 
splitters to aid in selection and blocking, and the possible 
use of a doubled- frequency , high-repetition-rate YAG laser. 

3C. Performance Predictions 

The theoretical performance of the LaRC Raman lidar can be 
estimated on the basis of equation (4) and the parameters in 
tables 4 and 5 assuming that interferences have been elimin- 
ated.” The questions of systematic error and calibration are dis- 
cussed below. Table 6 summarizes this predicted performance. 

When viewing 1000 ppm SO 2 at 300 m, the return signal should be 
0.55 counts per laser firing as registered by the DAQ in a 40- 
nsec bin (i.e., 6 m) . The accumulation of 100 counts should 
enable a measurement of SO 2 in the calibration tahk at night 
to a precision of 10% (i.e., a signal-to-noise ratio of 10). 

Such an accumulation would require about 200 laser firings or 
about 7 min. The potential for SO 2 monitoring at stacks is 
very promising. 

Diagnosis of correct operation and of freedom from inter- 
ferences of the Raman lidar involves the measurement of other 



Table 6, Performance Prediction for LaRC Raman Lidar for Various Gases. Systems 
parameters are listed in table 5 and include a 6-m range interval at 
300 m at a firing rate of 30 ppmin. 


Gas 

Concentration 

(ppm) 

Time Period for 
10% Measure 
(min) 

Return per 
Shot 
(cts) 

Remarks 

SO 2 

103 

7 

0.55 

— 

NO 

1Q3 

83 

0.04 

— 

CO 2 (7683) 

330 

76 

0.044 

— 

C02(7683) 

3 X 10*^ 

■ — 

— 

— 

O2 

2.1 X 105 

— 

1.23 

ND 1.3 added 

N 2 

7.9 X 105 

— 

0.72 

Other channel 
ND 1.0, etc . 



gases in addition to SO^. Table 6 contains entries for the 
predicted performance of the Raman lidar when viewing S02f 
NO, ambient CO 2 , high concentration CO 2 , ambient O 2 , and, in 
the reference channel, ambient N 2 . The N 2 counts are somewhat 
arbitrary depending on the configuration of the reference channel 
except that the 45 Mcps random count rate limit of the DAQ cannot 
be exceeded. Ambient O 2 operation is similar to the SO 2 opera- 
tion except for the change of the spectral filter and attenuation 
to keep the count rate below saturation. The configuration of 
each type of operation is outlined in conjunction with the dis- 
cussion of the channel transmissions and efficiencies appearing 
in table 7. 

3D. Review of Raman Lidars Used for Stack 
Plume Measurements 

A brief review is given here of all the Raman lidars used 
for stack plume or other pollutant gas measurements that have 
been described in the open literature as of the date of this 
report. Table 8 is a summary of this review. Assuming that 
the interference problem and lidar calibration have been 
successfully handled, the most important column for the pur- 
poses of comparing various lidar is the performance column. 

The measure of performance used here is the time period needed 
to make a measurement to a particular precision. This perform- 
ance time period is here determined by scaling the reported 
results of the various lidar groups to the hypothetical measure- 
ment of 1000 ppm SO 2 at a distance of 300 m with a precision of 
10 percent. A range interval of 6 m has been used and the use 
of a standard size 8-inch receiver has been assumed. The other 
colxamns of table 8 identify the Raman lidars and give several 
of their distinguishing features. 

As previously stated, Inaba and Kobayasi (ref. 7) first 
measured SO 2 in a nearby (30 m) oil smoke plume in 1970. Night 
measurements of SO 2 in a distant stack plume (200 m) were later 
reported by Nakahara et al. (ref. 14) . These workers were much 



Table 7. Transmissions and Efficiences of the Several Channels of the LaRC Raman 
Lidar. The lower half refers to the standard lamp class of operation 
and the upper to the Raman class of operation* 


Parameter 

SO 2 

CO 2 

O 2 

% 

^RAM 

0.025 

0.025 

0.0013 

0.00025 

^ratio 

0.91 

0.93 

0.94 

1.0 

RAM 

Tl 

1 

1 

0.05 

0.1 

^PMT 

0.16 

0.16 

0.16 

0.037 

^tele 

0.5 

0.5 

0.5 

0.5 


0.95 

0.95 

0.95 

0.7 (?) 

^atm 

0.9 

0.9 

0.9 

0.9 

■"sohott <5 

0.58 

0.61 

0.63 

0.46 (10 mm) 

filter 

0.66 

0.64 

0.69 

0.51 

X FWHM (A) 

27 

27 

28 

17 

-POl 

1! 

0.84/2 

0.84/2 

0.84/2 

0.88/2* ** 

sl 

Tl 

0.005 

0.005 

0.005 

0.01 

^sl 

5.8 X 10"5 

6.0 X 10" 5 

6.2 X 10 "5 

1.1 X 10"5 


* One way transmission over 300 m, 

** Depends on particular beamsplitter combination. Here leakage "adds” 1.05 factor. 


Table 8. 


Suiratiary and Comparison of the Performances of the Several Raman Lidars Used 
for Pollution Monitoring. Performance numbers are scaled from actual reports 
to a hypothetical measurement of 1000 ppm SO 2 at 300 m with a common 8-inch 
telescope. 


Remarks Reference 


Group 

Date 

Performance 

Inaba and Kobayasi 

1970 

400 min 

Nakahara et al. 

1972 

40 min 

Leonard 

1972 

33 min 

Hirschfeld et al. 

1973 

2.5 min 

DeLong 

1974 

1.5 min 

Melfi et al. 

1973 

6 min 

This paper 

1975 

15 min 

Kuper and Ebeling 


120 min 


6943 and 3371, spectrometer, 6,7 
plume at 30 m 

5320, filters, 14 

stack at 200 m 

3371, spectrometer, 20 

jet engine 

3472, spectrometer, 21 

daytime, gas clouds 

3472, etc., gas clouds 22 

6943, filter, 16,23 

stack at 200 m 

6943, filters/ 
cal. tank at 300 m 

6943, spectrometer. 


cal. tank at 50 m 
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more concerned with the demonstration of the technique and 
elimination of interferences than with performance although 
they did attempt an absolute lidar measurement. Their recom- 
mendations for the use of a high-repetition-rate , doubled- 

O 

frequency (e.g., 5320 A) YAG laser are somewhat altered today 
by the availability of efficient red photosurfaces and by NO 2 
fluorescence (ref. 18) and aerosol fluorescence (ref. 19) in 
the green. Interference filters were used to select the SO 2 
wavelength. The Japanese group also recognized the limitations 
of photon counting when a high background was present and when 
one had high instantaneous signal rates. Leonard (ref. 20) used 
Raman lidar to measure the composition and temperature of air- 
craft engine exhaust emissions. Among other features, he first 
used analog detection under computer control. Several descrip- 
tions of the Block Engineering lidar are available; see Schildkraut 
(ref. 10), Hirschfeld et al. (ref. 21), and Delong (ref. 22). 

This instrument is now based on a doubled-frequency ruby laser, 
a polychromator , and an analog detection system capable of day- 
time operation. It is probably the most powerful of all the 
Raman lidars due partly to its 36-inch receiver, but its complexity 
appears to limit its frequent use. Melfi et al. (ref. 16) used 
a Raman lidar based on the ruby laser and interference filters 
to measure SO 2 at a stack 200 m away. Melfi (ref. 23) indicates 
its larger telescope (i.e., 24-inch), but lower efficiency photo- 
detectors, would place this lidar about a factor of two better 
in performance than the present LaRC lidar. However, Melfi et al. 
(ref. 16) state that the quoted performance was obtained with an 
ND 0.9 attenuator in the return beam so that there is a discrep- 
ancy of about 2 to 4 between their results and the returns pre- 
dicted in table 6. It is possible that aerosol fluorescence 
might have caused this discrepancy (ref. 23) . 

In few of the above lidars was much attention paid to absolute 
performance. The reports contain very little information on which 
to base the performance estimates in table 8. 



None of the above groups measured a known concentration of pollu- 
tant gas in a calibration facility although many tried to estimate 
by other means the plume concentration. Measurements of nearby, 
artificial pollutant clouds were made, but mainly from the point 
of view of demonstration of the Raman technique. The present 
report is the first describing measurements of known pollutant 
concentrations and possible interferences under controlled 
conditions . 


3E. Possible Interferences 

3El. Introduction 

The design of the Raman lidar must include means to separate 
the various spectral components in the scattered light and mini- 
mize possible interferences. Two types of interference must be 
guarded against: leakage of the wrong wavelength return or 

acceptance of natural or induced broad-band background light 
into the spectral bandwidth of a receiver channel. The leakage 
can be due to either an inadequate rejection ratio of the spectral 
selector for strong returns at other wavelengths or to overlaps 
of the rotational structure of a Raman line with the spectral 
band selected. The broad-band background light could be either 
day background or fluorescence induced by the illuminating light 
pulse from gases or aerosols. If fluorescence is present, one 
may also have to correct for the additional absorption of light 
in the stack plume as indicated in equation (4) . If any of these 
items or the ones discussed in sections 3F and 3G are overlooked, 
the potential of Raman lidar summarized in table 5 cannot be 
achieved. 

3E2 . Filter Leakage Interference 

Filter leakage interference can be due to inadequate rejection 
ratios for Raman returns far from the line selected or for returns 
close to the line. Each are treated differently. Leakage inter- 
ference can also be due to an overlap of the rotational structure 
of a nearby strong line with the spectral band selected. One of 
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the most serious rejection problems is that of Mie scattered 
light, from plume aerosols at the same wavelength as the exciting 
line. Mie backscatter from atmospheric aerosols and diffuse 
scatter from the tank apertures and backstop are the corres- 
ponding problem at the calibration tank. An estimate of the 
seriousness of Mie backscatter can be made by using some of 
the typical backscatter functions. In the lidar geometry of 
table 5 , one would have an equivalent backscatter coefficient 
of 3 X 10“^° for Raman O 2 scatter, 2.4 x 10"® for Rayleigh 
scatter, and 1 x 10 ster"^ for a 10 km visual range ambient 
maritime aerosol. Assuming that the 3M black velvet paint on 
the backstop used in the tank tests is less than 0.5 percent 
reflecting, one can estimate its backscatter coefficient to be 
about 1.6 X 10"® ster"^. Backscatter from a stack plume might 
lie between that of the backstop and that from a white target 
with a value of about 0.3 ster"^ depending on the opacity of the 
plume. These estimates show that the Hie backscatter must be 
blocked from entering the spectral selection device by a very 
large factor. That factor might range from 10^^ to 10 The 
present Raman lidar makes use of interference filters (e.g., 
fig. 6) which are naturally suited to giving high transmission 

O 

for the 7546 A SO 2 line arid good rejection against Mie backscatter 

o o 

at 6943 A and N 2 Raman at 8283 A. The present filters were 
selected for relatively sharp shoulders and high transmission. 
Narrower filters can be inserted in the optical channel if found 
necessary. Other lidars make use of spectrometers for the selec- 
tion. These spectrometers have typically had excellent rejection 
ratios, but poorer transmissions. The use of spectrometers 
becomes attractive for lidars using lasers toward the UV where 
interference filters have poor transmissions and where lines are 

o o 

closely spaced (e.g.,, SO 2 at 3852 A and CO 2 at 3600 A) . The 
rejection ratios of both filters and spectrometers can be aug- 
mented by the use of absorbing glass filters. Figure 7 shows 
the transmission of a piece of Schott RGN-9 color glass used in 
the present lidar. A series of these glass filters with different 
thicknesses were obtained in order to optimize and study the 
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blocking of the Mie backscatter light at 6943 A. Tests of the 
efficacy of the blocking were carried out by observing the return 
from the backstop during the measurements at the calibration tank, 
and by observing returns from artificial plumes. Additional 
rejection capability without a serious transmission loss is 
probably best added to a filter Raman lidar by means of special 
beamsplitters used at low angles of incidence to obtain sharp 
shoulders in their transmission curves (e.g., fig. 8), One would 
transmit the weak return while reflecting the strong reference. 

The Raman lidar described here does not have such beamsplitters. 
Other strong returns distant from the selected line may break 
through the interference filter in addition to the strong Mie 

O 

backscatter. Figure 6 indicates that the N 2 return at 8283 A 
is also naturally blocked, but that blocking at longer wavelengths 

O 

may be a problem (e.g., water vapor at 9300 A). The latter block- 
ing can often be accomplished by taking advantage of the cut-off 
response of the photosurface used in the receiver. Figure 9 
indicates several Of these cut-off wavelengths. Use of a photo- 
surface yields much higher net sensitivity on-line than use of 
metal films to block out longer wavelengths. Note that the 
C31034 photomultiplier used in the present lidar is not sensitive 
to the water vapor line. 

Leakage of light through a filter due to a nearby Raman line 
(e,g, , O 2 or CO 2 for SO 2 work) depends on the design of the inter- 

O 

ference filter, Raman excitation at 6943 A is an advantage in 
that the Raman lines are well-spaced. This spacing allows one 
to use somewhat wider filters which have both good peak trans- 
missions and sharp shoulders. The width also minimizes the 
temperature dependence of the peak wavelength of the filter. 

The dependence of the peak wavelength on angle of incidence 
(ref. 25) is made use of in this lidar to peak the transmission 
and to verify that one is in fact viewing a Raman line rather 
than a broad-band return. Blocking of nearby returns depends 
on the sharpness of the shoulders of the filters. Figure 10 
shows a typical transmission curve of a three-cavity design. 
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all-dielectric filter which is currently available (ref. 26) . 

O 

For example, one can estimate the leakage through a 7683 A CO 2 

O 

filter due to the 7784 A O 2 Raman line. Figure 10 indicates 

O 

that a blocking of 2 x lO^* might be expected 100 A from the CO 2 
line center. Leakage at the 10 ppm level of CO 2 might then be 
expected. Section 5D presents measurements of ambient CO 2 which 
are relevant to this leakage question. On the basis of section 
5D, an SOp filter can be expected to block the O 2 Raman return 
to the 2 ppm level of SO 2 . Potentially more serious in an SO 2 

O <1 

measurement is the possibility of leakage due to the 7624 A CO 2 
line since the CO 2 in a pliime can reach 10 to 20% levels. Figure 
10 again indicates that leakage of Q-branch CO 2 scatter from 10 ^ 
ppm CO 2 in a plume is not serious (e.g., 20 ppm of SO 2 contri- 
bution) , 

Raman scatter, however, is not confined to a single line 
(i.e., the Q-branch) as is indicated in figure 11 for scatter 
from nitrogen. About 12% of the scattered light in the parallel 
polarization is distributed in branches due to rotational structure 
of the line. It might be possible for this weaker scatter closer 
to the filter pass band to contribute to the lidar return. If 
one assumes that the CO 2 structure is the same as N 2 for the sake 
of an order of magnitude estimate, one finds that convolution of 
figures 10 and 11 indicates only a 10 ppm contribution of 10 ^ ppm 
CO 2 to lO^ ppm SO 2 in a stack plume measurement. Similarly, one 
can estimate that the O 2 branches do not contribute to a CO 2 
signal at the 330 ppm CO 2 level. In view of the lack of explicit 
information about the CO 2 rotational structure and the uncertainty 
in the filter blocking for nearby wavelengths, it was decided to 
check for this type of leakage by viewing the calibration tank 
using the SO 2 filter when the CO 2 concentration was raised to 
10® ppm. Section 5G presents the results of this test. 

3E3. Broad-Band Background and Fluorescence Interference 

Two types of broad-band background can obscure the Raman sig- 
nal; natural background and induced fluorescence background. This 
interference can be minimized by proper spatial and spectral fil- 
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tering. Once minimized, it- can often be eliminated by suitable 
data recording and analysis techniques. Total elimination of 
day background is possible by working at night. Total elimina- 
tion of fluorescence is not possible unless a lidar wavelength 
can be selected that will not excite fluorescence. During tests 
with the Raman lidar at the calibration tank, it was found that 
day background rates were about ten times the standard lamp rates 
while looking at the black backstop. Measurements of N 2 and O 2 
Raman returns were found possible under these conditions if the 
day background remained constant and was subtracted by additional 
measurements. Operation of the sensitive SO 2 channel would not 
be possible, however, because the day background rate of about 
90 Mcps saturated the photon counting system. This saturation 

made data recording and subsequent background subtraction 
impossible. An analog data system would not be so limited (ref. 

21) . Induced fluorescence, if present, could also be subtracted 
from an observed signal. In this case, one would view the return 
at the selected line and at a line position nearby where no Raman 
return would be expected. A third channel could be used or a 
consecutive viewing period with the Raman line filter tilted 
so that it accepted only the diommy line. Again, one must minimize 
the fluorescence so that the photon counting system can record it. 
Induced fluorescence is known to be a problem at wavelengths 

O 

shorter than 6943 A with respect to both gases and aerosols (refs, 

O 

18, 19, 27). Little work has been done on this question at 6943 A. 
Melfi et al. (ref. 16) conclude that they observed no fluores- 
cence based on the fact that no decay time on the scale of 250 nsec 
was observed. They did not detune the filter, though, and a fluor- 
escence decay time of less than 250 nsec might have been present. 

Uo estimates of fluorescence interference are made here. Attempts 
to measure its importance are described in sections 51 and 5H in 
the context of the calibration tank. Definitive tests will have 
to I be carried out at a real stack plume. There one can employ 
at least three methods to identify fluorescence. The filters 
can be detuned to look for broad-band fluorescence, the range 
intervals can be narrowed to identify decay time constants as 
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short as 60 nsec, and, in the case of aerosols, the stack pre- 
cipitators can be used to vary the amount of aerosols independently 
of the gases. 

3F. Tests for Correct Operation 

Tests for correct operation of the lidar typically consist 
of viewing a standard lamp placed at the range interval of interest 
and measuring N 2 and Q 2 profiles through the atmosphere. These 
tests insure that the lidar data system is working correctly and 
that all optical elements of the lidar are in good condition. 
Correct operation of the data system means that the photomul- 
tiplier plateau has been reached, that the photomultiplier gate 
is properly adjusted, that multiple pulsing at the discriminator 
is not present, that noise does not enter the system during laser 
firing, and that the count rate saturation values are known and 
worked within. The first three questions as well as other photo- 
multiplier related questions are discussed in Appendix 2. The 
latter two questions are discussed in sections 5A and 5B. 

The optical performance of the Raman lidar can be checked by 
measuring an N 2 or O 2 profile in the atmosphere. However, it is 
just as effective and somewhat easier to view a standard lamp 
placed down range. Furthermore, the standard lamp can be used 
to diagnose any efficiency loss that might be noted in the N 2 
tests. The predictions of Raman lidar performance in table 5 
were not entirely achieved. The search for the discrepancy 
using the standard lamp is discussed in section 5E. The lamp 
used was a 1000 W GE DXW Quartzline lamp with a nominal spectral 

O 

flux density at 7500 A of 23 liW/cm^ nm at 50 cm (ref, 28). The 
photon flux falling on the 8-inch aperture telescope at 300 m 
is approximately 9 x 10^ photons per sec-A and is denoted by B. 

O 

The photon flux at 8283 A is greater by 1,04. The count rates 
expected from the standard lamp in each channel are given by 
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where it is assumed that B and the system efficiency e are 
constant over the small spectral width of the filter. The nor- 
malized spectral width of the filter is given by 
for estimation purposes is just the full-width half maximum of 

O 

the filter in A. The final factor, x, represents any attenu- 
ation added which is not normally used. The system efficiency 
can be expressed as 


^s£ "^atm '^tele '^BS '^Schott '^pol ”^p ^PMT ”^1 


( 6 ) 


where 


^atm i^ep^ssents the transmissivity of the atmosphere from 


the lamp, '’•tele represents the transmission of the telescope, 

Tnc represents the transmission (or reflectivity) of the beam 
splitter for the polarization in use, "^poi represents the trans- 
mission of the polarizer, ''^schott represents the transmission 
of the blocking filters, x^ represents the peak transmission 
of the spectral filter, represents the quantum efficiency 

of the photomultiplier, and xj represents the transmission of 
any other optical attenuator normally used in the optical train. 
All of these transmissions and efficiencies were measured at 


either the wavelength of operation or one close by. Filter 
transmissions and shapes were measured on a Cary 17 spectro- 
photometer. Photomultiplier efficiencies were measured using the 
photomultiplier as a diode, an HeNe laser, a laser power meter, 
and the relative sensitivity curves of the photosurface. All 
of these parameters were determined in conjunction with the 
absolute efficiency search described in section 5E and are 
listed in table 8. The parameters in this table are arranged so 
that the components of the efficiency for Raman operation can be 
found by reading up from the lower dashed line and the components 
of the efficiency for standard lamp operation can be found by 
reading down from the upper dashed line. Table 9 lists the 
expected standard lamp rates for the parameters of table 8. 

Note the slight correction for polarizer leakage in standard 
lamp operation. Section 5E compares these predictions with 
actual measurements and discusses any discrepancies. 
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Table 9. Predicted Rates in the Several Channels of the 
Raman Lidar Viewing the Standard Lamp. 


Channel 

Rate (Mcps) 

SO 2 

14 

CO 2 

15 

O 2 

16 

N 2 

1.9 


The standard lamp also aids in the optical alignment of the 
lidar. Using it, one can position the collimating lens and the 
lens that focuses the light on the small photosurface of the 
RCA C31034. One Can also check for vignetting by moving the lamp 
through the field of view of the telescope. Auxiliary spectral 
lamps placed close to the telescope were also used. These were 
an independent check on the angle-tuning constants of each filter 
and on the location of the peak transmission wavelength. 

3G. Calibration 

The difficulty of making an absolute lidar measurement makes 
it imperative that a suitable calibration method be developed. 
Many workers have suggested using a second reference channel, 
usually N 2 , to aid in this calibration. The concentration of 
the unknown gas can then be written in terms of ratios of known 
or measurable quantities. For example, the ratio of N 2 to O 2 
concentration can be written as 



where the wavelength ratios are given in table 8, the cross section 
ratios in table 4, the efficiency ratio is determined by the standard 
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lamp measurement, and the transmissivity ratio accounts for any 
optical attenuation added to or removed from individual channels 
for Raman operation if different than for standard lamp operation 
The final factor is the ratio of observed range interval counts 
corrected if need be for pileup. It is here assumed that the 
filters are peaked for each Raman line. It is obvious that all 
the ratios in equation ( 7 ) are known except for the relative 
efficiency of the two channels. Similarly, any other gas con- 
centration can be determined relative to N2 using equation ( 7 ) 
as long as relevant relative efficiences can be determined. 

The determination of the relative sensitivity of the two 
channels can be done in several ways. Measurement of the N2/O2 
ratio in the atmosphere (i.e., 3 . 7 ) could serve to measure the 
efficiency ratio at these two wavelengths in addition to deter- 
mining atmospheric and plume transmissivities and to pi ividing 
laser pulse energy normalization. If the Raman line of the 
unknown concentration gas is close to the O2 line, no further 
measurements would be necessary for calibration. In fact, one 
may wish to measure the ratio of unknown to O2 as reference in 
this particular case. If the measurements are consecutive in 
time with a single channel, no linking of channels would be 
necessary. If done simultaneously in a dual channel mode, the 
two channels must be linked using Raman O2 measurements in both. 
If the Raman line of the unknown gas is far from the O2 line, 
the spectral gap must be bridged by use of the standard lamp. 
Either of the two reference methods outlined above could be used. 
Work in this report emphasized the use of the standard lamp and 
the N2/O2 ratio calibration methods. Most of the results quoted 
in section 5 were obtained with standard lamp calibrations, A 
number of these results were recalculated using the N2/O2 method 
as quoted in section 51 . 

The relative Raman efficiences can be obtained from the 
relative standard lamp efficiencies using equation ( 5 ) 



and using equation (6) 

^RAM ^sZ ”^atmy^^pol 

The most convenient form of the relative Raman efficiency 
formula is 

O2 

^RAM 

N2 

^RAM 



where the peak transmission of the filters allow the use of a, 
the ratio of filter areas as determined by weighing. If the polar- 
izer is used for Raman measurements, only the 1.05 factor remains 
for N 2 reference measurements. Otherwise, the full factor must 
be used for N 2 reference measurements. O 2 reference work would 
not have this correction factor. The subscripts of equation (10) 
must be changed to correspond to the actual case at hand. The 
two basic calibration formulas are equation (7) and equation (10). 

The relative calibrations just discussed yield absolute con- 
centrations of pollutants to a precision determined by the length 
of time available to fire the laser. A loss of efficiency in a 
channel will, not affect the ability to make the measurement, but 
only the time it will take to obtain enough returns to get the 
desired precision. The optimization of system efficiency is 
discussed in sections 3F and 5E. The accuracy of the measure- 
ment depends on the calibration procedure. A measure of the 
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lidar accuracy can be obtained by comparing lidar measurements 
of known concentration pollutants in the calibration tank with 
in situ measurements. These me^asurements are reported in section 5. 

4. CALIBRATION TANK 
4A. General 

A Raman lidar calibration facility has been constructed .to 
contain the gases for the calibration measurements. The tank 
consists basically of a steel tube 2 m in diameter and 20 m 
long. In principle of operation, the tank is charged with a 
quantity of gas which is mixed with ambient air to a known 
concentration. Raman lidar measurements can then be made 
through a known volume and concentration of gas. 

4B. Tank Description 

The calibration tank consists of a 2-m’-diameter by 20-m- 
long steel tube. The tank is charged with gas through a 15-m- 
long perforated manifold. Two fans, located in the tank, are 
used to mix the charge. At each end of the gas tank there is 
an aperture 0.6 m in diameter. The aperture is required to 
maintain the gas mixture in the tank. A larger diameter aper- 
ture or no aperture would greatly reduce the gas residence time 
in the tank. Each aperture is fitted with a remotely operated 
shutter. The shutters are closed during short periods when 
measurements are not being made, thus increasing the gas resi- 
dence time. A 1.3-m-wide air curtain fan is located at each end 
of the tank. Any outside air currents blowing down the tank 
would result in a loss of the gas charge in the tank. The heating 
and cooling system is used to regulate temperature and for dehumidi- 
fication. This capability has been used primarily, to date, to 
dry the system for use with SO 2 . Figure 12 is a schematic of 
the tank showing its principal components. 
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4C. Gas Operation 


The tank has been operated for Raman lidar calibrations with 
various concentrations of CO 2 and SO 2 . Gas mixtures have been 
maintained to within ± 5 percent of a desired level for residence 
times of 10 minutes. The tank has also been charged with varying 
concentrations of CO without difficulty. 

4D, Aerosol Operation 

Aerosol screens have been produced in front of the tank to 
simulate a smoke stack plume. The aerosol screens were produced 
by burning a mixture of used automotive engine oil and varsol. 

A waste cotton material Was used as a wick-base for the smudge 
pot technique. The smoke was then vented through the air curtain 
fan to produce a screen of aerosols about 1.3 m wide by 15 cm 
thick. 

4E. Instrumentation 

The calibration tank instrumentation (table 10) consists 
primarily of gas analysis and. temperature measurement. The hand- 
ling of the gas sample for analysis is made with an Environmetrics 
gas handling pump. The gas sample is taken about mid-way in the 
tank through 6-mm-diameter teflon tubing. The sample is then 
pumped through the analyzer at a rate of about 1 to 3 standard 
cubic feet per hour. A chart recorder is used to record the 
gas concentration with time. Figure 13 shows a typical time 
history of the gas concentration for a CO 2 and an SO 2 run. 

Table 10. Gas Sampling Instrumentation 

Gas Manufacturer 

CO MSI LIRA 

CO 2 MSI LIRA 

SO 2 Environmetrics 


Range 

0 to 1000 ppm 
0 to 20% 

0 to 5000 ppm 
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5. LIDAR MEASUREMENTS OF AMBIENT GASES AND HIGH CONCENTRATION 
POLLUTANTS IN THE CALIBRATION TANK 

5A. N 2 /O 2 Ratios in the Atmosphere 

Measurements of the N 2 /O 2 concentration ratios serve to 
check for correct operation of the lidar, calibrate the lidar 
for; absolute measurements of pollutants, and to measure the 
atmospheric extinction. Initial measurements with the LaRC 
Raman lidar during these tests were concerned with the check 
for correct operation. The latter two possibilities were pursued 
during measurements through the tank. Figure 14 is a display of 
an early return from the atmosphere alone. The counts in each 
bin of the sixteen-bin photon counting system accumulated after 
a number of laser-firings were displayed here after the range 
square correction of equation (4) was made. The error bars on 
each point are statistical precisions. The N 2 and O 2 profiles 
were obtained consecutively since the sixteen bins form only a 
single channel. Range bins of 200 nsec were used for a total 
range time of 3.55 ysec or 511 m. The returns for each gas should 
show an exponential decay depending on the amount of extinction 
present. The scales are arbitrary in that the returns were atten- 
uated to make them about equal and to place the initial count rates 

within the limits of the photon counting system. Figure 14 indi- 
cates* that the fields of view of receiver and transmitter do not 
overlap until about 1.55 ysec for the receiver field of view 

used for this run. Not so obvious is the need for a pileup 

correction at the beginning of the profile to account for the 
random count limitations of the photon counting system (e.g,, 

45 Mops) . Figure 15 shows the ratio of the above profiles. The 
Units are again arbitrary along the abscissa since the relative 
sensitivity of the two channels is not known without additional 
information. Note that the adjustment of both returns to about 
the same value in figure 14 corrects for much of the overlap and 
pileup problems without explicit calculation. The ratio curve 
in figure 15 should be a straight line with range. Appendix 1 
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compiles the data analysis programs that allow presentations of 
data and results such as in figure 14 and figure 15 along with 
the various corrections (e.g,, pileup) . 

5B. N 2 /O 2 Ratios Through the Calibration Tank 

Most of the Raman lidar work was done in conjunction with the 
calibration tank. Initial N 2 /O 2 profiles and ratios were taken to 
verify that backscatter from the tank and tank apertures would not 
break through the blocking into the Raman channels. Elimination 
of this backscatter was found to be necessary for the highest 
sensitivity channel. It was accomplished by proper configuration 
and blackening of the tank and by proper choice of auxiliary 
blocking filters. Figure 16 shows the atmospheric profiles 
through the tank on the same range scale as figure 14. The 
profile is truncated by the backstop that was placed to make 
the lidar range safe to passers-by. Note that the nitrogen 
channel shows no sign of the backstop whereas the oxygen channel 
shows its presence. The oxygen channel is obviously poor in 
blocking the Mie scatter from the backstop even though a 
5-mm Schott filter was used. The interference filters in this 
test were soon replaced with new filters. Neither channel indi* 
cates the presence of the tank which normally appears about 1.9 
psec. It was found that apertures as small as two feet in diameter 
could be used to help contain gas in the tank as long as the laser 
beam was collimated to 1 raillirad in full width, the tank apertures 
were kept smooth, and kept black with 3M Black Velvet paint. 

Tests with the highest sensitivity channel are discussed in 
sections 5D and 5E. The residual counts after the backstop 
in figure 16 are due to afterpulsing in the photomultiplier 
which is discussed in Appendix 2. Both tubes were RCA 8852 
at this time. 

Once it was proven that the Raman lidar could look through 
the tank successfully, attention was placed on the atmospheric 
profiles again for purposes of exploring the calibration of the 
lidar and measuring atmospheric extinction. Figure 17 displays 
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sample returns from the atmosphere through the tank. The time 
bins are now 40 nsec which allows two bins to be placed within 
the tank as discussed in section 5c. Figure 17 shows no evidence 
of tank breakthrough. The backstop is just out of range. Figure 
18 shows the N 2 /O 2 count ratio accompanied by an absolute scale 
determined by standard lamp measurements. The observed ratio of 
3.8 ± 0.1 indicates fair agreement within the expected 3.7. It 
is typical of N 2 /O 2 profiles. This calibration proceeded as dis- 
cussed in section 3G with the properties of the interference 
filters known by laboratory measurements and with their peaks 
adjusted to the Raman line wavelengths empirically by tilting. 
Standard lamp ratios were first taken with an auxiliary counter 
connected to the photon counting system. This method was later 
replaced by taking the calibration data with the complete system 
itself (triggered by an auxiliary pulser) and by analyzing it 
in a manner similar to the profile ratios as discussed in Appendix 
1, Such a standard lamp profile is a good check on correct oper- 
ation of the detector. The fluctuation of the data points in each 
bin of figure 17 raises questions concerning the true width of the 
time bins, the true randomness of the counts, and the possibility 
of noise during laser firing. The latter is the most difficult 
question to answer due to the low Raman return. It was eliminated 
by observing the standard lamp through a SO 2 filter with no SO 2 
gas present and with the laser firing out of the receiver field 
of view. The first two questions are closely intertwined. Exten- 
sive tests and analysis with the standard lamp showed that fluctu- 
ations such as appear in figure 17 are nomal and that the widths 
of the bins were correct to a few percent. 

5C. High Concentration CO 2 

The encouraging results with the Raman lidar for atmospheric 
constituents prompted measurements of other gases in the cali- 
bration tank. CO 2 was chosen because it is not toxic like SO 2 . 

A range of concentrations was chosen that spanned a return signal 
as small as that from 1000 ppm SO 2 up to ten times that signal. 
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The CO 2 concentrations ranged from 10,000 to 35,000 ppm. Figure 19 
shows the CO 2 profile and the N 2 reference profile. The CO 2 return 
was taken with the more sensitive RCA 31034, The channels were 
taken consecutively in time. The tirae scale of operation (e.g., 

100 shots) does not allow the ambient 330 ppm CO 2 to be seen. 
Operation of the CO 2 channel without CO 2 in the tank is a good 
means to monitor the condition of tlaa black surfaces of the tank 
and any resultant increase in Mie backscatter breakthrough from 
the tank. Routine use of this test is time-consuming, however, 
so it was replaced by an analog Mie detector installed as a third 
channel in the Raman lidar detector package. Tank Mie returns 
were typically 0.01 times the backstop return. To obtain absolute 
CO 2 concentrations from the lidar measurements shown in figure 19, 
the standard lamp calibration method was again used. Figure 20 
compares the Raman lidar measurements of the various CO 2 concen- 
trations with the concentrations measured by the in situ gas . 
sampler. The systematic error appears to be smaller than about 
10 percent. Similar results were obtained during several other 
high concentration CO 2 tests. The Raman lidar could thus be 
used to characterize the high CO 2 concentrations in a stack plume 
if those were of interest (ref. 23) . These results with CO 2 indi- 
cated that the Raman lidar, the calibration tank, and the lidar 
calibration method were ready for the more demanding measurements 
of 500 ppm SO 2 . The absolute return of the 19 December 1974 test 
corresponded to 0,7 counts/shot which is a factor of 3.4 lower than 
predicted by table 6 once the presence of a second 5-mm Schott 
filter in the CO 2 was taken into account. The standard lamp CO 2 
rate was 5.8 MHz which is a factor of 1.6 smaller than predicted 
by table 9 once the extra Schott filter is taken into account. 

These discrepancies in absolute return and performance instituted 
a lost photon hunt which is summarized below in section 5E, 

5D. Ambient CO 2 

Table 6 indicates that a Raman lidar capable of detection and 
measurement of ambient CO 2 at 330 ppm should certainly be capable 
of measuring 5C0 to 1000 ppm of SO 2 as long as interferences are 
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not present. The time period for such a measurement would be 
twelve times that for the SO 2 measurement or a predicted 76 minutes. 
In view of this long period and the capability of placing known 
concentrations of gases in the calibration tank, no work was 
done specifically on ambient CO 2 . However, in the course of 
the high concentration CO 2 measurements, signals from ambient 
CO 2 were obtained from those time bins outside the tank. These 
signals imply that the ambient CO 2 concentration was 850 (± 170) 
ppm which is much too high unless leakage of CO 2 from the tank 
raised the concentration to that value temporarily. Assuming 
that gas leakage from the tank could not cause such a high reading, 
one can set an upper limit for the blocking of the CO 2 filter to 
©2 Raman scatter, Mie backscatter is not expected since two 5-mm 
Schott filters were in use. The excess counts when attributed 
to ©2 Raman scatter indicate that the C ©2 filter blocking for 
the ©2 line (100 A away) is only a factor of 1.4 x 10^. Figure 10 
indicated that one might have expected a blocking of 2 x 10*+. 
Narrower filters could be ordered to obtain a higher blocking 
if desired. Since the CO 2 , © 2 # and S ©2 filters are all very 
similar in peak wavelength and shape, one can estimate the 
blocking expected for these filters and other lines as discussed 
in section 3E. ©ne should try two C ©2 filters in tandem to try 
to eliminate the leakage, 

0 

5E. Absolute Efficiency Search 

The high concentration CO 2 results of December 1974 indicated 
that the standard lamp signal was low by a factor of 1.6 and the 
C ©2 Raman signal low by an additional factor of 2.1 from the pre- 
dicted signals based on the parameters of table 7. The situation 
was not clear at the time because not all the system parameters 
were then as well known as implied in the discussion of table 7. 

The 3,4 discrepancy was serious since it made the difference 
between 6 min and 20 min in predicted performance. The search 
also took longer than might be expected because of modification 
of the single channel photon counting system to a dual channel 
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system, installation of gating provisions for the C31034 photo- 
multiplier, studies of partial reduction in auxiliary blocking 
in the O 2 channel, and a variety of other changes and mishaps. 
During this time period, the extensive statistical tests were 
carried out. The particular C31034 in use from December 1974 
to May 1975 probably degraded in use by a factor of two over 
this period and was, in fact, destroyed by accident on 2 May. 

A concurrent problem v;ith the photomultiplier bias circuit also 
took place on 2 May. After installation of a new base and tube, 
the LaRC Raman lidar performed at its highest sensitivity. A 
N 2 /O 2 profile on 16 May gave returns of 0.45 counts/shot for N 2 
and 0.60 counts/shot for O 2 . The standard lamp rates were 9.6 
Mcps for O 2 and 1.4 Mcps for N 2 . Comparisons to table 6 and table 
9 indicate that the standard lamp rates are down by a factor of 
about 1.5 and the Raman rates are down by a factor of 1.3. The 
airfield visibility that night was reported as 6 miles in haze, 
clearing after a rain shower to 7 miles. The light rain that 
night did not appear to affect the Raman returns. The slope of 
the N 2 and O 2 Raman returns showed little evidence of any absorp- 
tion due to the atmospheric conditions. One concludes that the 
discrepancies are real. Extensive tests of the receiver were 
carried out during which field of view, focusing, aberrations, 
photosurface scans, etc. were all studied. It is suggested that 
the receiver discrepancy of 1.5 may still be due to a poor tele- 
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scope mirror reflectivity in the 7600 A region in addition to the 

O 

50 percent transmission measured at 6328 A. If this reflectivity 
is the sole cause of the receiver loss, the Raman return can be 
increased by a factor of 4 by recoating the mirrors. The dis- 
crepancy in the Raman rates by a factor of 1.3 cannot now be 
explained. It may be due just to uncertainty in one of the 
other transmissions or efficiences in table ?. 

5F. SO 2 Measurements 

During May to July, four observations of SO 2 concentrations 
in the calibration tank were made. Only two of these were made 
after peak performance of the lidar had been reached. Figure 21 
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shows the results of the analysis of the measurements of 19 May 
1975. The range profile is not shovm since it has all the same 
features as the CO 2 range profile in figure 19. Most of the 
data for 19 May was taken in the dual channel mode so that a 
simultaneous N 2 normalization profile was obtained. The cost of 
this could have been a systematic error due to the single SO 2 
range’ interval not being centered in the tank. The numbers at 
each point in figure 21 are the number of laser firings needed 
to obtain that point with its corresponding precision. This 
result in addition to the other SO 2 runs is the basis for the 
conclusions of this paper that the LaRC Raman lidar is capable 
of night measurements to 10 percent accuracy of 1000 ppm SO 2 in 
stack plumes at a distance of 300 m if interferences are eliminated 
or accounted for. 

In these tests, the SO 2 was certainly being detected because 
the signal except for a small residual disappeared when the gas 
was released from the tank. In addition, the SO 2 filter was de- 
tuned by tilting and the signal also disappeared. 

5G. Gas Interferences 

Several possible gas interferences with the measurements of 
the Raman lidar have been discussed in section 3E. One of these, 

O 2 Raman leakage through the CO 2 filter, was examined in section 
5D above. That examination indicated that the inteference filters 
appeared to perform 20 times worse than the expected behavior 
given in figure 10, and led to estimates in section 3E of other 
leakage possibilities. For example, the O 2 Raman return should 
be blocked to the 2 ppm level when using the SO 2 filter. The 
stronger N 2 Raman signal should contribute even less. The closer 
CO 2 Raman line at 7624 A was expected to contribute about 20 ppm 
to an SO 2 measurement if the CO 2 concentration was 10 5 ppm. Pre- 
liminary results of this CO 2 test of the SO 2 filter on 26 June 
1975 was that the leakage signal could not be distinguished from 
leakage through the SO 2 filter in the absence of CO 2 in the tank. 
The background breakthrough of the SO 2 filter corresponded to 
0.01 counts/shot or 40 ppm of SO 2 . The cause of this leakage is 


42 



not known. The only other prevalent atmospheric constituent is 
water vapor, but the photomultiplier photosurface is not sensitive 
to its wavelength. One should try two SO 2 filters in tandem to 
help distinguish the cause of this leakage. 


No tests of gas fluorescence were made other than the detuning 
of the SO 2 filter when SO 2 was in the tank. One may wish to test 
NO 2 for fluorescence leakage (since this is one noise source in 
the green, tef. 18) in spite of the fact that its concentration 
will be quite small right at the stack (tens of ppm) . 


5H. Particulate Interferences 

Mie backscatter interferences and possible particle fluores- 
cence interferences can affect the Raman lidar results as dis- 
cussed in section 3E. Mie backscatter from stack aerosols is 
the possible problem at the stack. Mie backscatter from atmos- 
pheric aerosols and diffuse scatter from the tank apertures and 
backstop are the problem at the tank. The residual counts in 
the SO 2 Raman channel might be a problem of Mie backscatter from 
cunbient aerosols. It will be of interest to check this hypothesis 
here by. scaling from, the observed backstop breakthrough signal 
using the estimated backscatter functions of section 3E3. Similar 
means to check and eliminate Mie backscatter would be used at a 
stack, A stack plume is not expected to be as an effective 
scatterer of light as a black (or white) target, especially with 
the present emission regulations* The observed return from the 
black backstop using the SO 2 filter and a 6-mm RGN-9 absorber 
is about the equivalent of 3 x 10^ ppm SO 2 . The tank is not 
usually seen under these conditions. Assuming that the black 
backstop is still less than 0.5% reflecting (ref. 29), one can 
estimate a return of about 150 ppm SO 2 for a typical aerosol 
breakthrough. A residual SO 2 signal of about 30 ppm SO 2 was 
observed. However, the analog channel showed that scatter from 
the teink apertures was of greater magnitude than aerosol scatter 
at the range of the tank so that the tank apertures should also 
have been seen in the residual SO 2 signal. Also, addition of an 
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extra 1 iran of RGN-9 did not eliminate the residual SO 2 return as 
it should have done. 

Fluorescence of aerosols is a problem that can only be studied 
at a stack due to the difficulty of imitating the aerosols in a 
real plume. The authors found it difficult to produce a plume 
just for the purpose of checking the blocking that would be needed 
at a stack much less one for the study of fluorescence. At least 
another 5-mm RGN-9 Schott glass filter will be needed at a stack. 

It is possible that the residual SO 2 signal is due to ambient 
aerosol fluorescence. 

51. Alternate Calibration Schemes and Systematic Errors 

The alternate calibration scheme to the use of a nitrogen 
reference channel is the use of an oxygen reference channel. Its 
advantage is its proximity in wavelength to the SO 2 Raman line. 

The spectral lamp need not span such a wide response region. In 
addition, the identical channel can be used with only the use of 
an auxiliary optical attenuator to keep the count rate comparable. 
Until the lidar is modified to have two identical channels, the 
method requires consecutive profiles. A test of this scheme 
directly has not yet been carried out due to other priorities 
and equipment difficulties. However, a number of both CO 2 /N 2 
a;nd SO 2 /N 2 measurements were made the same night as O 2 /N 2 ratios. 

It was then possible to determine the SO 2 concentration, for 
example, in two ways. The nitrogen normalization technique plus 
standard lamp yielded 1130 ± 140 ppm. Use of the N 2 /O 2 ratio in 
addition to the N 2 /SO 2 ratio to obtain an absolute SO 2 concen- 
tration yielded 1150 ppm SO 2 . The data run was number 7 of 19 
May. The sampling indicated 1367 ± 257 ppm. The systematic 
error was not affected in this hybrid calibration technique. 

Whether or not the direct O 2 calibration would eliminate the 
systematic error is not obvious. A similar test with CO 2 data 
of 10 December 1974 indicated an increase from 76,000 ppm to 
84,000 ppm using the O 2 hybrid scheme. The actual in situ CO 2 
concentration was 74,000 ppm in this case. Most of the N 2 /O 2 
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ratios were also high, again indicating a systematic error and a 
lower than expected O 2 concentration. The cause of this system- 
atic error is not known. 

6. LIDAR MEASUREMENTS OF AEROSOL PLUMES AND ATMOSPHERIC 
EXTINCTION AT THE CALIBRATION TANK 


6A. Theory and Review of Other Lidar Measurements 

Many workers have suggested the use of remote Raman scatter 
to obtain the transmittance of the atmosphere and aerosol plumes. 
Leonard and Caputo (ref, 30) have made the most thorough investi- 
gation of such possibilities experimentally for the atmosphere 
and artificial smoke clouds. In the present context, one would 
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view the range dependence of the nitrogen Raman scatter at 8283 A. 
Equation (4) indicates that this range dependence would lead to 
a measure of transmittance since the cross section and distri- 
bution of nitrogen is known. Figure 22 shows schematically such 
a return for many laser firings. The slope of the line depends 
on the atmospheric transmittance as in figure 17 and the discon- 
tinuity on the plume transmittance. The smoke plume transmitt- 
ance can thus be determined in principle as long as single scatter 
prevails and enough radiation traverses the plume. It is also 
assigned that there are no other absorbers present and that inter- 
ferences such as fluorescences have been eliminated or accounted 
for. The principle of this transmittance determination using the 
Raman lidar is. very similar to the lidar method tried by Cook 
et al. (ref. 31) using Mie scatter from ambient aerosol before 
and after the plume. The Raman scatter method has several 
potential advantages which will alleviate difficulties with the 
Mie scatter method. 

For many purposes, knowledge of the atmospheric or plume 
transmission is sufficient. For example, a model for radiative 
transfer through the atmosphere including aerosols may only require 
the net transmittance to be measured. For plumes, one of the 
Federal Emission Standards (ref, 2) limits the opacity of the 
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plume to 20 percent or less except for 2-min periods every hour. 

The opacity measurement is also important for the proper interp- 
retation of SO 2 measurements in the plume, Leonard and Caputo 
(ref. 30) measured cloud transmittance to several percent pre- 
cision over a range of one-quarter mile. The potential for the 
LaRC Raman lidar for plume transmission measurements will be 
discussed below. 

The second Federal Standard (ref, 2) for particulates in stack 
emissions requires the mass emission rate to be below 0.18 g/Mcalj^. 
Reference to table 1 indicates that s percent coal, for example, 
will emit 8.2 s g/Mcalj^ input. A plant burning one percent fly 
ash coal would therefore need a precipitator operating at 97.8 
percent efficiency to enable the plant to meet the emission 
standards. Again, the standard is not written in terms of rates 
so that all plants must have the same efficiency precipitators 
and the larger plants can put a larger absolute amount of par- 
ticulates into the air. The Federal Standards for emissions 
(ref. 2) specify the methods to be used to determine the mass 
emission rate of particulates. Due to difficulties of location, 
operation, etc. of these other devices, it is useful to raise 
the question of whether or not a transmissivity measurement could 
be interpreted in terms of a particulate mass emission rate. 

Conner (ref. 32) has reviewed this question as well as the 
various methods for measuring opacity and transmissivity of 
plumes. There is evidence in certain cases that opacity can be 
rather directly related to particulate mass concentration in a 
plume. This mass concentration must be related to mass emission 
rate by means of supplementary measurements of plume exit velocity 
and diameter. 

Mie lidar measurements of stack pliimes yield another parameter 
which depends on the aerosols and their concentrations. That 
parameter is the Mie backscatter coefficient, ' which would 

replace the na in equation (4) . The lidar return, however, now 
depends on both 3j^(tot) and these cannot be deter- 

mined exactly unless they can be otherwise related. Collis and 


46 



Uthe (ref, 33) address this problem and show how 3,(7 t) can be 

** 

approximated based on additional assumptions and/or measurements. 
One such measurement is just the transmission one already dis- 
cussed. Another is to calculate the relationship based on 
particle size data obtained by sampling, Johnson and Uthe (ref. 
34) obtained measures of plume (length-wise) mass content for a 
stationary source stack in this manner. A third method would be 
to find an empirical relationship between the backscatter coeffi- 
cient and the opacity. Collis and Uthe (ref. 33) obtained encour- 
aging results in a controlled scattering chamber with injected 
fly ash. They found indications that a lidar wavelength of 0.7 
ym yielded results dependent on opacity whereas a wavelength of 
1.06 ym yielded results more closely related to mass concentra- 
tion. It is not clear that any of these aerosol backscatter 
measurements possess any advantage over the opacity measurement 
discussed in the prior paragraphs. 

6B. Predicted Performance of the LaRC Raman Lidar 

The LaRC Raman lidar can be used either in a single channel 
mode or in the dual channel mode to measure plume transmittance 
by recording the nitrogen signal as a function of range. Figure 
22 displays schematically such a return. Note that the back- 
ground atmospheric transmission is also measured as the slope 
of the general curve. The plume transmittance is measured by 
the discontinuity at the plume location. For estimation pur- 
poses, one can reduce the analysis to the data points from the 
range bin before and after the stack. The transmission of the 
plvime can be written as 

N 2 /N 1 = T = e”''^ 

where N 2 and Ni are the counts after and before the plume 
and T is the two-way optical thickness of the plume. Assuming 
the absence or correction for systematic errors, one can express 
the error in the optical thickness as 
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At = / [e^ + 1 


on the basis of random statistics of the received counts. The 
error in transmission can be expressed as 

AT = TAt 

Table 11 is a summary of the time periods expected for the stated 
number of total counts needed in the first bin to obtain 5 percent 
measures of high T and 10 percent measures of low T. It has 
been assumed that the N 2 channel attenuation has been adjusted to 
yield about 1 count/laser firing. Higher returns are handled 
with increasing difficulty by the photon counting DAQ. Night 
measurements are also assumed. Note that the time periods get 
long at both high transmissions and low transmissions. No improve 
ment can be made with the LaRC lidar unless a high-repetition-rate 
laser is used with less attenuation (or a larger receiver) or 
unless an analog DAQ is installed with the same proviso. The 
latter case would allow day operation with proper choice of 
telescope size. 

Possible interferences are fluorescence and inadequate 
blocking, both of which have been adequately treated in section 
5. The advantage over the conventional lidar method of measuring 
transmission (e.g.. Cook et al. , ref. 31) is that the intense 
return from the plume is blocked before reaching the detector 
and does not cause overloading of the analog DAQ or afterpulsing 
of the photomultiplier. Both of these problems are serious for 
a conventional Mie lidar, 

6C. Results 

Several attempts to measure pliame transmittance of fake 
plumes at the calibration tank during the interference studies 
were not successful due to the difficulty of producing the type 
of plume needed for the interference studies. It is recommended 
that the plume measurements be continued during a field trip to 
a stack. 


Table 11. Time Periods Needed to Measure Various Plume Transmittance Using the LaRC 
Raman Lidar to View Nitrogen Raman Scatter. Returns in that channel have 
been adjusted to yield 1 count/ shot in the interval before the plume. 


Transmis s ion 
T 

Optical Depth 

T 

Ni 

(counts) 

Ax 

AT 

Period 

(min) 

0.90 

0.1 

1300 

0.04 

0.04 

43 

0.82 

0.2 

900 

0.05 

0.04 

30 

0.37 

1.0 

400 

0.1 

0.04 

13 

0.13 

2.0 

600 

0.1 

0.013 

20 


Monitoring of atmospheric transmission was quite difficult 
due to the short baselines involved for both Raman and standard 
lamp measurements and due to the low aerosol attenuations. 

Over the period 18 December 1974 to 10 July 1975 many nitrogen 
profiles were obtained. They showed a general correlation with 
the highest lamp rates (e.g. , 1.7 Mcps) corresponding to flat 
range-squared nitrogen profiles and with the lowest lamp rates 
(e.g., 1.52 Mcps) corresponding to Mie attenuation coefficients 
of 0.5 Km“^. The standard lamp extinction over 300 m is not 
inconsistent with such a coefficient (e.g., 0.86). The two-way 
Raman extinction would be 0.75 for the worst case. The typical 
Mie coefficient for the nights of operation was 0.25 Km"^ or 
0.86 two-way Raman extinction. The extinction of the nitrogen 
profile in figure 17 was also the greater value, but did not 
fit into the above sequence due to the use of a different lamp. 
Kildal and Byer (ref. 12) quote a Mie extinction coefficient of 
0.24 Km"^ for visibilities of 10 km maritime haze. Visibilities 
in the visible obtained from the nearby airfield were 7 to 10 
miles on several of these typical nights. These visual visibili- 
ties were determined by observers in the flight tower as they 
looked at nearby features. Both observer and feature were probably 
above the ground haze which was common on the marshy lidar range 
in early evening. None of the measures of Raman return or stand- 
ard lamp rates were taken for the express purpose of determining 
atmospheric extinction. They are therefore not accurate enough 
for; determining the low extinctions that were typical of the 
nights of operation. Closer attention to the measurement would 
greatly decrease the errors and provide a more accurate measure 
of atmospheric extinction. 

7. CONCLUSIONS AND RECOMMENDATIONS 


The tests of the LaRC Raman lidar at the calibration tank 
at a distance of 300 m indicate that night measurements of SO 2 
concentrations in stack plumes are possible. Accuracies of 10 
percent are achievable within 30 min integration times for 500 ppm 
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of SO 2 at the stack exit. Comparable measurements of HCl at a 
rubbish-fueled stack should also be possible. Supplementary 
measurements are required to obtain the concentration and mass 
emission rate. The possible interference of aerosol fluorescence 
can only be evaluated at a real stack. If fluorescence is present, 
it can be subtracted from the signal as long as reasonable signal- 
to-noise ratios exist. Simultaneous measurements of plume trans- 
mission are possible with this dual-channel lidar. Accuracies of 
better than 10 percent in transmission of the aerosol plumes should 
be possible within 30 min with the LaRC lidar. An early test at 
an instrumented stack is strongly recommended. 

The time period for the above measurements can be decreased 
somewhat by improving the telescope transmission and finding the 
source of the absolute discrepancy with predictions. With the 
present lidar and photon counting DAQ, these time periods can be 
decreased by about a factor of 4 until the bin count rates approach 
random count saturation. A larger telescope or larger laser pulse 
energy would not help once this saturation rate is exceeded. For 
night operation, a high repetition-rate laser at the same or higher 
average power in conjunction with a larger telescope would decrease 
the time periods of observation even more. The high repetition 
rate Nd:YAG laser is known to excite the fluorescence, however. 

To make day operation possible, one really needs to have the com- 
bination of a somewhat larger telescope and a laser pulse of 
higher energy. The present lidar with its efficiency discrepancies 
eliminated should be capable of day operation if an analog DAQ 
were added in the several channels. The analog DAQ should possess 
the capability of the same 20 to 40 nsec time-bin operation. 

Analog operation would also aid the plume transmittance measure- 
ments since these rates were artificially lowered in the present 
lidar and aid the discrimination against fluorescence. If the 
new DAQ were added, one may also want to consider the use of an 
O 2 reference channel rather than an N 2 channel and make both 
channels identical. The size of the lidar van prevented this 
step in the present lidar. The advantages of the O 2 reference 
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are that one does not need a special beamsplitter (a piece of 
glass will do) and the spectral distance between the SO2 line 
and the O2 line is quite small, making standard lamp calibration 
more reliable. In HCl operation* one would still use the N2 
reference channel. 

In conclusion, the present lidar is an adequate instrument 
for monitoring SO2 and plume opacity from stationary sources at 
night. Its large depth of field is a definite plus. Its compact 
size and ease of maintenance and operation compare favorably with 
larger, spectrometer-based Raman lidars* The N2 channel maintained 
its performance over the whole year of operation. Its photon 
counting system makes possible this small size and provides very 
modest but adequate data-recording capabilities. That photon 
counting system also restricts the lidar to night operation. An 
analog DAQ added to the present lidar in addition to the improve- 
ment of receiver efficiency should provide a powerful device for 
the remote measurement of SO2 and HCl stationary source emissions 
at the 500 ppm level and for the measurement of plume transmissi-> 
vities during day and night. 
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Figure 1, Vibrational Raman Scattering Leonard, 1972). 


cn 

cn 



i 



AX = X^ Al'/C 


O- (7T) »( 1/q - 

i 

Figure 2, Raman Shifts in Frequency and .Wavelength. 
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Figure 3. Raman Backscatter Function (Theoretical) for Typical Oil Smoke Plume 
Using Ruby Lidar (after Inaba and Kobayasi, 1972) . SO 2 Concentration 
is 1000 ppm. Positions of Scatter from HCl and Liquid Water are Also 
Indicated. 
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Figure 4. Laser and Receiver of NASA-Langley Raman Lidar. 




Figure 5. Data Acquisition System of NASA- Langley Lidar. 





Figure 6. Typical Performanc 
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Figure 10. Typical Behavio^ of Three-Period Interference 
Filter with 30 A Full Width at Half Maximum, 
(Ditric Optics) . ■ 
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Figure 13. Calibration .Tank Performance for Several Gases. 
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Figure 15. Ratio of the Range Corrected N 2 and O 2 Profiles in the Atmosphere. 
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Figure 17. Raman Signal from N 2 and O 2 Through the Calibration 
Tank with 40 nsec Bins. 
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Figure 18. Ratio of the Range Corrected N 2 and O 2 Profiles Through 
the Calibration Tank. 
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Figure 19. Raman Signal from High Concentration CO 2 in the Calibration 
Tank. Nitrogen Reference Signal is Included. 
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SYSTEM PARAMETERS 
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Figure 20. Comparison of CO 2 Raman Lidar Measurements with In Situ 
Gas Analysis. 
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Figure 22. Schematic of Nitrogen Raman Return Expected When Viewing 
Through a Stack Plume. 
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DATA ACQUISITION OPTIONS AND DATA ANALYSIS PROGRAMS 
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MEMORANDUM 


SUBJECT: 

Raman LIDAR Data Reduction Program 

DATE: 

1 August 1975 

TO: 

Dr. 6. B. Northam, NASA/ESSD 

FILE: 

V-19000/5NAS- 208 

FROM: 

Ms. C. F. Jones, LTV/HTC 

■ 


CC: 

Mr. M. L. Brumfield, NASA/ESSD 
Mr. T. E. Chappell, LTV/HTC 
Dr. S. Poultney, ODU 




I. INTRODUCTION 

Data from the Raman LIDAR project can be taken in single or dual channel 
modes. Two classes of LIDAR operations occur: a pulsed-laser class which 
is used to monitor Raman backscatter from one or more types of gases; and 
a standard lamp class which is used to calibrate the LIDAR. The data 
* reduction program was designed to be able to handle any combination of these 
modes and classes by keying on several input parameters. In the single 
channel mode, there are numbers in 16 bins corresponding to data return as 
a function of time. If a single mode is used, that mo^ is usually repeated 
for and paired with another possible channel of operation. The numbers in 
each bin are usually the result of an accumulation of counts over a certain 
time period of operation (here measured in laser firings). The program will 
average compatible sets of such operation for up to two gas types and adjust 
for the range-squared dependence in the LIDAR equation-^fov' the laser-pulsed 
class of operation. Error bars are calculated for each point. The resulting 
data is fit with an exponential whose exponent is a measure of atmospheric 
extinction in the case of laser operation on an. atmospheric constituent and 
then plotted. With paired single channel modes or with dual channel modes, 
this data is then ratioed, plotted on a separate graph with error bars, and 
fit with a horizontal line. All standard lamp measurements must, be analyzed 
omitting the range-squared correction, of course. In the dual mode, the 
16 bin single channel is split into two interleaved parts; with different gas 
types or standard lamp wavelengths recorded at the same laser firing. The 
program proceeds as above for each of the interleaved parts except that points 
for the ratio must be picked off from the best fit exponentials prior to 
‘ ratioing. 
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II. INPUT 


The first data card feeds general information to the program. An 
example appears in Appendix I. The input constants are as follows: 


BT - Initial time delay (usee) between laser firing 
and data acquisition. 

TAU - Time width (usee) of successive bins 

F - A pile-up correction factor 

TYPEl and TYPE2 - Types of gas to be analyzed 

IDUAL - Code = 1 for single channel mode; = 2 for dual channel. 

ISTL - Code = 1 for laser firing mode; = 1 for standard lamp 

mode. 

DATE - Month/day/year on which the data was taken. 

The following data cards contain the actual counts as well as the total 
number of laser firings used to produce those counts. Different gas types 
are input on separate data cards. If operation is under dual channel mode, 
blanks are left on each data card for those alternate bins which do not apply 
to that particular gas. .Also included on these data cards are data, run 
number and gas type for easy identification purposes. 
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III. 


DATA REDUCTION 


A. Pile-up Correction 

Data cards are read in and count rates from each bin are corrected 
for system pi 1 e-up error: 

^bs = ^eal 

where = observed count rate (I.e., counts/, 

width) 

Freai “ actual count rate after pile-up correction 
f = rate correction term (presently equals 126 MHz) 




A 


The pile-up correction subroutine (TRANS) uses the Newton-Raphson 
iterative method to solve the above equation within an absolute error 
of .001. After a solution is found, the actual rates are converted 
back to counts and transfer is made back to the main program. Corrected 
bin counts for each gas type are then accumulated into two arrays (CNTl 
and CNT2). 


B. Averaging and Range Correcting 

The last card in the data deck is a blank card which signals the program 
that there is no more data. Data is then averaged if needed due to multiple 
runs. If in laser pulsed mode, range corrected: 


R. = BT + (t/2) + (i-1) XT i = 1, 16 

where R. = range of each of 16 bins (usee) 

BT = begin time 
T = time width 


R.^CTN, = R? (Z CNT 
^ ’ j=l ij 


N 

/Z NOS ) i 
j=l j 


= 1, 16 


where CNT. . 

* J 

NOS. 

U 

N 


= counts of a particular gas (i ranges through the bins 
and j ranges through the different runs) 

= total number of shots for run 
= total number of runs 


i 
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And the Poisson error for a particular gas type and bin is: 


n n n 

V £ CNT n NOS 
h ’ ’ j=i ji j=i j 


Program checks the standard lamp parameter ISTL, and if turned on 
sets R? to one to cancel the range correction factor in the above 
and ^ following equations. 


C. Exponential Fit 

Exponentials are then fitted to the data using a least squares 
technique. Blank bins (either due to system overflow or to dual 
channel operation) are not included in the fit. If we define: 

In R^CNT. 
m p p 

{Z E i=l,m m < 16 

then a library subprogram will fit a linear curve by minimizing the 
sum of the residuals: 

m ‘ m o 

I = EW [(a + a X ) - Y ] . 
i=li i=li 0 1 i i 



Once the exponentials have been calculated, blank values are replaced 
by the fit values. Standard error for the exponential fits are determined 
by the sum of the residuals: 



m 

/(Z 

i=li 


) (m-1) 



D. Ratioing 

Ratios of TYPEl to TYPE2 gases and deviations on the rHios are 
calculated: 

R12. = CNT, / crifp 
* i ^i 

^R12. = VVCNT^, + I/CFT 2 R12^ i=l, 16 
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Finally, the weighted ratio mean and its standard error are given by: 



IV. OUTPUT 

Computer print-out includes data input, pile-up corrected data, range 
corrected and averaged data and ratio data. Plot output includes curves 
for both types of gases as well as their exponential fits. If two gases 
are to be analyzed a ratio plot is also drawn. Also printed out on the 
plots are the exponents and ratio with their estimates of error. 




83 


■ DATA CARO SPECIFICATIOMS 


Card Ho. 1 

Variable Name 

BT 

TAU 

F 

TYPEl 
TYPE2 
I DUAL 

ISTL 

DATE 


Description . Columns 

Begin time (usee) 1-6 

Time width (usee) 7-12 

Pile-up correction (126 MHz) 13-18 

First type gas (e.g., N2, 02, C02, etc) 19-22 

Second type gas (leave blank if none) 23 - 26 

Code =1 for single channel mode 

. = 2 for dual channel mode 28 

Code = 0 or blank for laser pulsed mode 

= 1 for standard lamp mode 30 

Date of data (mpnth/day/year) 33 - 40 


(Note: For variables TYPEl and TYPE2, if less than 4 letters, must be 
right justified. For example, N2 would be in columns 21 and 22 if TYPEl.) 


EXAMPLES 


Single Channel Laser Pulsed Mode (IDUAL=1, 1STL=0) 
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Single Channel Standard Lamp Mode (IDUAL=1, ISTL=1) i 
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Dual Channel Standard Lamp Mode (IDUAL=2, ISTL=1) 
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Cards Mo. 2 



Variable Names 

Description 

Columns 

DAYS 

Date of data {month-day-year no 
slashes 

1 - 6 

RUN 

Run number (right justified) 

7 - 8 

TYPE 

Type gas (must match either TYPEl or 
TYPE2, right justified) 

9-12 

NOS 

Number of shots for that run (right 
justified) 

13 - 16 

DATA(i), i=l, 16 

Counts from the 16 channels (each 
channel is alloted 4 columns on data 
card, right justified) 

17 - 80 

NOTE: Any number of this type data card is allowed. Program will 

average 


any number of runs of a particular gas. Different gas types must be 
on separate cards. Leave blank columns for overflow data and for 
alternate channels on dual channel data. Last data card must be 
followed by a blank card. 










Sample data card set-up 



3 na aoo ?a 


Gas type 2 data 
cards 


f 030 ?? 6 l a Q 3 aoo 76 


/•■jaOrrBl 1 ns SCO 95 


f 030 ?^|:.l 3 - Ha aoo 


(oaorrSi a Ha aoo 


H ) 30 ?.'PI 1 HS aoo 


1 . SjO 

“ 1 1 — T . 


. 0-; 1S6. M3 D3 3 0 03/0?/?5 


; FORTRAN S rAT.EMENT I UJCHTiriCAVIUti 

»iC»l*C* * ' ^ I 

‘oluTTO 0 liSlIGCDOaDDOSOIiaOCIIDOODDOIlllDlI gTToI [IOQOQIIDODOSOSO!!ODCD!!9000I)QDOOOOOO|8 0 0 0 0 0 0 0 
I 1 t 1 I ; • t II II 1/ II II n It II II I) II I! II I) :i ii n ii ini n ii ii :i ii n it ii n ii •: ii non o <( •' ii « » ii :i » » is » si » » u ii t: n ii » it ii u is ii ii ii ii it is it i) is n ti 

Ml 1 1 1 11 1 1 1 n n t ] t n 1 1 1 1 11 1 1 n It I M 1 1 1 1 1 ! 11 11 1 n n n 11 II M 1 1 n 1 1 n n It n 1 1 n 11 n 


2<2 2 2 212 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 7 . 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 


M 2 3 2 n 3 0 3 ) 2 n 3 2 3 3 3 3 3 2 3 0 3 2 3 3 1 3 3 3 3 3 3 3 3 0 3 3 3 3 3 3 3 3 3 3 3 3 2 2 3 3 } 3 3 2 2 2 3 ) 3 3 2 1 3 3 2 3 3 3 3 2 0 3 3 2 3 2 

I ■ ' 

4>U 4 ( 4 4 < 4 4 4 « 4 4 4 4 4 4 U 4 4 3 4 4 4 < 4 4 4 4 4 U 0 -{ 4 4 4 4 4 < 4 4 4 4 4 4 44 4 4 4 4 4 4 4 M 4 4 < 4 ( 4 « t 4 U 4 4 4 4 4 4 4 4 4 4 

i’O 2 3 5 3 S 3 $ 3 3 S $ S S S 3 S S 3 $ 3 S S 5 S S 3 S S 3 S 3 ‘ 3 9 3 5 5 S S S S S 3 S 9 2 9 5 S S S 9 S S 0 S : S S S S 3 S S 5 S S S 3 S 3 0 S S S 2 S 0 

I 

(is t 6 6-’0 0 0 S 0 C S S S C 6 0 $ S S 6 6 0 0 0 6 G 0 C e 0 8 6 S 0 0 0 6 6 t S t 3 3 C S 6 0 0 0 0 C S 5 0 S 0 0 E 0 0 0 9 S B G 0 6 0 0 S 0 1 S { 0 9 3 S 0 

3 ' 

7I1 T 7 ) ; I ? n 7 1 77 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 n 7 7 7 7 7 7 7 n 7 1 7 7 7 n 7 7 T 7 7 7 ; 7 77 7 7 7 7 7 n 77 7 7 7 7 7 7 7 77 7 7 n 7 

1 

' (ii 3 3 c 0 1 8 n ( 1 1 e 1 8 0 0 s 8 8 « 3 1 8 B 0 0 CO c I e I s ( e t ( 1 1 ( G 8 e n 8 { 0 1 ( 1 8 a t e 1 1 0 8 0 s ( c 1 0 ( ! t > 8 0 0 B a 3 ( 8 1 

' ' , ' ’ 

' o'a 9 3 9 oio S 3 8 3 8 9 9 9 8 9 9 9 3 9 8 9 8 8 9 8 3 9 o 3 9 8 9 3 8 9 9 9 8 3 9 9 ! 9 9 S 9 9 3 0 8 9 0 3 8 8 8 9 9 0 3 9 9 3 8 8 n 3 J 0 0 0 3 8 3 8 9 8 

>, I 1 !’ 1 I :|t;' ! I til! III! II IS It 11 11 ii i!ii II o II IS isii It II H II 11 DM »» !• msuii III1II 0 II H It n SI II SI nun it II iisi is ii ii tiu ts it n tin ii ii it ii » i»it n ii iii: 

9 * Ml«y 








Rflnos 


SAMPLE OUTPUT 


ORTE 05/02/75 

(Note: Dual path range refers to the propagation to Mnnc 1 

the hackitop end return.) riL/Ub. 1 



■'DUflL PRTH RANGE (MICROSECONDS) 


• X lO-i 
ffiP* 



N2 

EXPlfil + 02 X RANGE) 
01 =' . 1392 : 

02 = -.1363 

STONDflRD ERRORr ,0302 


02 

EXPffil + R 2 X Rft/vJfiEl 
fll= .0532 

fi 2 = -.0912 . 

STRNDRRD ERRORr ‘.0162 





. DUAL PATH RANGE ( M ICROSECGNOS ) 


87 





ORTE 06/25/75 
MODE 2 

STD. LRMP CPLIBRRTIOM 


X 10-“ 



X 10*‘ 
61“ 


OJ 


EXPCfil + fl2 X RRNGei 
Rl= -S-itOS? 

R2= - . ! >i87 

STRNDRRD ERROR= .0108 


S—S- 


03 


-e- 


■ET 


4h- 


m 


-a- 


m 


-a — 3 — B- 


3 


-e — ^r-as 


02 • ' 

EXPCni + R2 X RRNGE) 
Rl= -.e^iso 
R2= ' .OOOii 

STR\0RRD £RRC:fe .OilS 


lp~ . 


- 0 - e — 0 — @ — © — © — Q- 


gf — &- 


-te*; 


X IC- 


cu 

'■»o 


-2 


OL 


I ! 


i _j_ 


1.50 1.5: 


I."”? 1.75 1.80 1.85 1.80 1,35 2. S 2.05 2.10 2-15 2.&;i 2.'.?5 ; 2 .-n 2,35 2.4,0 2.4: 

{^iCR0.8EC.a\'2SJ 


89 



VO 

o 




. RAMAN LIDAR DATA REDUCTION PROGRAM 

• 



NcGRAV PCULT(IAPLffCLTPTjT) 

* 

LCCCC2 


C 1 M F.N S I ON ■ 7 1 Pi 4 ) , S I G ( A ) 


C C v; G C 3 


i: IPENSiriN YUH, n,X( lEl ,RIrS 1 C (in, 1 }| SUMSOn )tA( 2t 2) ,R(Z| 1) tC( 2|1) 


CCCGC'J 


' UIiXfcNSION CNTK III} , CN 72 ( 1 H ) , R ( 1 H ) , EG 1 { 18 ) , EG2 1 10 ) , R1 2 ( 1 8 J , CATA ( 3 A ) 


, f FR12 (10 ) tWK ua,2 )tVJl IB) |V<12 no) 

CCCCC2 


IiNTEGFfl BUN 



C 

INITIALIZE VARIAOLES TO ZirRO 


cccccb” 


NCSl=NOS2-C . 


CCC0C5 


i:g 1 1 = 1, u 


CCCCC7 


CMl I n=CNT2 ( I ) = 5U2 (I )=EGI { n=EG2( n = ERl2( n = G. 


cccc n 

1 

CC1N7INUF; ‘ . 


(u:qv>2c 


REAC 100, HT,7AL,F,7YPF.l,TYPE2,inUAL. ISTL, CATE 


cccc -'i ■•r~ 

' ICC 

rCs'.XA7 { 3F£. 2 ,2AA ,212 ,2X, A» ) 


GCO.i-|/, 


■ PiilNT ■ 202 , iiT ,'7AU,'F , 7 YPF 1 , T YPF2 


C C 0 1 ' t- 2 

2C2 

FCRPAt i 1HI,=:‘CATA I NPUT«/ IX , <■ IN IT I AL CEL AY=4 , FS .3/ IX , >>P.GSOLUTrCN = « 

• 



, ,Fd. 3/lX,*F (PILL UPl =«,FQ.371 X,*GA's TYPES >:',A1C,* ANC*,Al6) 


CCCCfc2 


FT l'^'r=BT + (TAU/2 , ) 


CC0Cc2 

?. 

REAO 101, DAYS, FUNjTYFe ,NCS, ICATA ( n , 1=1 ,16 ) 


'COOIC? 

ICl 

FCRNAT { A6 , 12, AA, lA, 16FA.C) 


CCC1C7 


IF (P-LN.tO.OJGCK 7 


CC'JllC 


' i’lilM r 203', DATf; , RUN', T YP E , Nil S , ( DATA"! 11,1 = 1,16) 


CC0122 

2C3 

FCPMAT ( IX ,'1'CAT E’:* ,A0 , « fUJN<',I2,’> TYPE*,A6,=!' NG . SHDTS^ , I A , AX , 1 6 F 




,i;-G) 


CCC 122 


1 F ( TYPF.FC.TYPCI 13,5 



C ACCCf^L-LAT5 TYPEl GAS 


GC0127 


ACS 1 = NQS 1 + NO.S 


CCC 1/. 1 


r:c A 1=1,16 

• 

rcoi‘<i2 


IF (DATA! n.EQ.C. 1 GGTO A 


CCCIA3 


r-i:.3s=DATA ( n 


:jcc 


CALL TRANS 1 FQO S , FR HA L , F , TAL ,NCS ) 


cccioc 


CATA(n = FREAL 


CCC1C2 


C.NTII n = CKTl( I ) +FRL-AL 


C'JCU-i 

A 

CCMINUF 


f;c;c;iC7 


PRINT 212, (HATAl n ,1=1,16) ' 


Ci;oi7i 


(-GiU 2 


CcCi )'i 
t C C C 1 1 s 

C ACCOLLATI; TYPC2 GAS „ 

5 NGS2 = KriS2fKQS . . 

'..xc,A....i=i.r.it '1 : 

IF (lUTAi t ).cQ.c.) c.nr 

% 

' i;CtH IL 


i-rfis=L'ATA ( n 


, (;cc2cc 

. . CCfj2C2 

. 

. ...G^U.aRAN!) .(.i\Cfl!i..r.f.Rl;.AJ..t.F.aAU.,AC^5.1.. , 

. „ .LAU (. U = FJA.EAL 







C CC 3 i ^ I F_ij_S T L . 1.C . l_L i£J AI_2 0 7jlT_Y_P EJ 

ecu 3 2 " ~ P>. INT 26 5, (C N’THI ) , I = L»_i V 

C C C3j 4 ~ FEINT 20 ^, (E G I ( I ) t 1 = 1 . 16 ) 

0 ~C(r 2 ^.C ;CA ECRmat 1 lHC«5«;y,lCH p EUAfiCN£/- l;Li-Ld£^^ 

C C C 3 5 C 2C 5 _ FIJLfM T 1.LLQ i.tXj. tJLe.*_i 

j 7X,*0o,7X,»H»»7X.»I»t7 Xt»J*t 7X .■> K» . 7 X . «L ■ > ,7 X,* > . 7 X . t 7 X . »n<« , 7 X » 


CCTr'SC 

----- 

- 

'pfKf'A f tVr/)32Xt*«ANCE CilRREXLTGIL_^VjF_^JS CJ±xA6J . - 

0C03 50 

1C 7 

FCRMATCAC/ n2X.'«STAMnARO LAMP AVERAGF'? HF^.Atl _ _ . - 

q_C C15 r 


IF (ISTL.EC.O) FPIM 206.TYPE2 

CC0_15 7 


IF (ISTL.EC.l) PRINT 2C7.TYPF2 - - - 

CCC26 7 


FPIM 205, (CM2( !)♦ 1 = 1. 16) . ... 

CCC^C 1 


FPl.VT 20 A. ( FG2 ( ! ) .1 = 1 .16 ) __ _ . 

CCOA 12 


FFRwaT (20X.#AFTFR P II tlJP CH PR GC T I CN * , 7 X , t 6 F ‘ - 11 


f. FIT 

EXPfNFNTIAL 7C MTHCGcn cATA 

CCUi- 13 


i-T<;*i=o. 

CCOA 1 A 



C C 0 A 1 5 


LT 120 1 1 = 1.16 - - 

CCOi It 

1201 

W.FLW=v^?;UR♦FG) n »*’>? -- - - - - -- 

CCCA22 


1 = C . 

ccca:-2 


l'c 1202 1 = 116 . 

f CCA?5 


IF I CMl ( I 1 - re. C. ) GfTr 1202 

• C C 0 A J t 



”ICCCA’C„. 


Hi J = 1\5 u G LQ.J ♦ *1 .. 


VO 



CCOA32 


Y iL ) = ALCG(CNT1( in 


CCCA24 


>(L)=Pm 


CCC^^C 


VTSV = k>(L )*V»TSF 

• 

CCCA-i2 

1202 

CCNTINUE 




ITCT = L 

* 

CC0/»A6 


CALL LSOPCL tl«,LTOT,X, l,Y,W,2,2tRFSICfSUFSC,A,n,V«K,IERP) 


C FILL I\ ZFKCcS IN C^TA wiTH EXPGFENTIAL FIT (TYFfl) 

CCCAd3 


CC 12C 1=1,16 


CCC^£5 


IF ICNTl 1 I l.NF.C. ) GCTC 120 


CC04£6 


CNTl ( I ) = EXP( R( 1 ) +B( 2 ) *K ( I ) ) 


CCC^?5 


TSC = R ( I ) ♦*2 


CCU't7 7 


IF (ISTL.EC.l) TSL=l. 


CC05C2 


tP12( I1 = ER12( imSQ/lCNTK n*NCSl) 


CCC5C7 

12C 

-12( I)=CNT1( I ) 


CC0513 


SIJVSC=S(JRT (SUPSC/(WTS.v*(LTCT-l ) ) ) 


CCC523 


CALL FSELCC 


CCCJ2^ 


CALL CALPLT tl.,1.,-3) 


C POIM EXPHNEMIAL LATA CF FLCT 

CCC527 


CALL NTTATE ( 1 2. , 6. , . 2 , TYPE 1 ,C . ) 


CC0533 


CALL NCTATE ( 1 2 . 5 , 5 . 6 , . 1 5 , 20F c XP ( A 1 ♦ A2 > RA^GE ) , C. , 2C ) 


CCC527 


CALL NCTATE ( 12 . 5 , 5 . 2 , . 1 5 , 5H A 1= ,0.,5) 


CC05<i3 


CALL NUMBER ( 1 3 , 5 , 5 . 2 , . 1 5 , 0 ( 1 ) , 0, , 4 ) 


CCC5A7 


CALL NCTATE ( 1 2 . 3 , A . 8 , . 15 , 5F A 2= ,0-,5J 


0C05 53 


CALL NUMBER I 13. 5, A. E,. 15, 01 21,0., A) 


CC0557 


CALL NCTATE ( 1 2 .5 , A . A , . 15 , 1 5 F S T AND A PC EP RCR= , C . , 1 5 ) 


CCC5£2 


CALL NUMBER ( 1 A . 5 , A . A , . 15 , S LF S C ,0 . .A ) 


ccns£7 


IF (NCS2.EC.0J GCTC 322 



C FIT 

EXPONENTIAL TC CXYGEN DATA 


CC057C 


V'TSM = SUMS0=O. ^ 


CCC572 


► SL!^=0. OW 


CC0572 


cn 12 C 6 1 = 1.16 


CCC574 


CESICII)=0. 


CCCSTf 

12C6 

V-5LM = FSUM*-EG2n )*«2 9,? 


CCOftCl 


L=0 9) |p4 


CCC< C2 


nc 1207 1*1,16 ip ^ 

• 

CCCdC4 


IF (CNT2I I l.EQ.C. 1 GCTC 1207 'd ^ 


LCCL C5 


L = L + l Kp 


CCCdC7 


v(L)=WSUM/EG2( I )*^2 


CCOd,l2 


Y{Li = ALCG(CNT 2 n n ZATo 


CCCd 15 


ML)=R(I) 


C cot 17 


FT<;p = Fn. )4 -wt';m 


CCCfc2 l_ 

-U0 7 . 

CCNTINUE . _ 


.CCCC.2.3 


■j-Tn.T=L.. . _ -- - 



i 


ccctzt 

CALL LSCPCL (18.LTCT,X,l,Y,V<,2.2.PESrC,SUf'SC.A.C.KK. lERP J 

. 


C FILL IN d"RnCS IN CATA WITH fXPnNFNTIAL FIT tTYPF 2 J 


CCCfe«2 

rC 122 1=1.16 


ccot^-^ 

IF (CNT2( I l.NE.C. ) GCTC 121 

. 

CCCd^S 

CNT2 ( I) =EXP( C(I ) *0(2 ) *R( I n 

• 

CC065'< 

T£g=RU )**2 


■■ CCOtid 

IF (ISTL.Et.l) TS«=l. 


;' ccceu 

5P12m=ER12( I)+TSC/(CN'T2(I )«NCS2) 


CCO'btt 

121 F 12 ( H = R 12(1 l/CNT 2 m 


ccctn 

FP12( M^Sf3RT{tR12{ I ) )*R12 m 

• 

CCOd77 

122 continue 


CCC7C1 

SLNSC = SQRT (SUNSC/tHTSP^ILTCT-l ) 1 1 


CCO'» 1C 

PRINT 20 e 


CC0713 

2Ce FCOVAT (4(/)32X,*RATlCS OF TYFEl TC TYPE2 GAS*) 


'CCC711 

PRINT 205 , IW12( I) .1 = 1.16) 

• 

ccc7:s 

PRINT 20^, IEM12( 1). 1 = 1.16) 


C PPIN7 C-XPMNFNTI/U ffllA CN PICT _ .. 

003727 

call NOTATE 1 1 2 . . 3 . . . 2. T Y PE 2 . C . . ^ ) 


CC07A2 

call NCTATE (l2.£,2.d..l5,2CFEXP A1 4 A2 X R AN 0^= ) . C . , 20 ) 


' CC07A7 

CALL NOTATE I 1 2 . 5 , 2 . 2 , . 1 5 . 5H A 1 = .0.,5) 


CCC 7^2 

CALL NUMBER U J . 5 . 2 . 2 . . 15, C ( 1 ) . 0 . , A ) 

-r ; 

CCC 757 

call notate ( 12.5 , 1.0 15,5HA2» .C.,5) 



C'J07£3 

CALL NUMBER 1 1 2 . 5 , 1 . F , . 1 5 . C ! 2 ) . 0 . . A ) 


CCC7£7 

CALL NCTATF ( 1 2 . 5 , 1 . A , . 15 , 1 5 F S T A NC A P C ERROR = , 0 15 ) 


CC0773 

CALL NUMBER ( 1 A . 5 . 1 . A . . 15 , SUR SC . 0 . . A ) 


CCC7 77 

22 ? call ASCALE (R. 10. .16. 1.10.) 


CUCC2 

IIMI 3)=XCKG = RI 17) 


CCICC5 

^ I«(A ) = XS»0( IB) 



CCIC I 1 

CALL AXES(C-.n-. 0 .. 10 ..XCRn .XS.-E.C.fPFCFANNP! S..15.-P1 




C PUT iLL TYPEl ANC TYPE2 GASES IN CNE ARRAY ANC CALCULATE ORIGIN. ANC 

■ILALE-.- 

CClCiA 

J = 0 


coic;5 

r.C IICO 1 = 1.16 


CCICP7 

CATAIJ41 )»CNTl( I ) 


001021 

IF INCS2.EC.C) CNT2I I )»CNT1( I) 


C0103«i 

CATAIJ + Slx CNT2 ( I ) 


ccjc;*7 

J = J42 


r, c 1 c ^ 0 

HOC CONTINUE 


CCIC^ 2 ., 

c c i C « d 

CALL ASCALF (OAT A ,6 . , 1 . 1 0 . ) 


CNTK 1 7)=C\T2( 17) = SIG(3) = YCPG=DATA (3 2) 


CCU5m 

CNT 1 ( IS ) =CNT2( IC )=S IC lA | = YSC = 0 AT A ( 3 A ) 


Q£:i 'i' 1 

CAll AXP5 I 0. , C. . 90 . , 6 - , YOPC ,Y*50 . ) - .0. .TYPE) . . 15. A ) 

, 

CC1C7^. 

call axes ll0.jC.,.90,,6...yCRG.YS!l, l-jC..iJ_YPE2,.l£..-A.l 



_j;l ! iLni GAS typel 




CCIIC7 

123 CALL LINPLT{R,CNT1, 16,1. -1.11.1.0) 



001 1 n 

cc no 1 * 1,16 



CCl 121 

I 1 ) = T I«(2 ) = R ( I J 



CCllJA 

nGm=cNTi( i)-EGi(i ) 



CCllic 

S IG(2 ) = CNTH I )+ £GU I J 



CC113C 

CALL LINPLT (TIF,SIG,2,l,0,0,.l,0) 



CCllAC 

lie CCMTINUF 



C PLCT EXPCNC\*TI/»L FIT TC Q6S TYPE 1 

CCl 142 

nc 12C5 1=1,16 



OCl l<i<i 

12C5 CNTK n = EXP{B(l )+Hl2)<Rl I » ) 



CCllf 5 

CALL LINPLT ( R , C M 1 , 1 6 , 1 , 0 , 1 1 . 1 ,0 ) 

C ' 

!r_ 


UCl 16^ 

IF (N0S2.EC.0) GCTO GS 




C PLCT GAS TYPE2 

6 g 


C0116S 

CALL LINPLT ( R , CN T2 , 1 6 , 1 , - 1 , 1 2 , 1 ,0 ) 

Ss 


CCl 175 

CC 12 1=1,16 



CCl 177 

TIV(1)=TIM2)=RII) 



CC12C2 

SIG(1) = CNT2( I1-£G2{I ) 



CC12CA 

£IG(2)=CNT2( I ) + EG2lI ) 



CC12C6 

12 CALL LINPLT ( T I P , S IG , 2 , 1 , 0 , C , . 1 . C ) 




C PLCT EXPCNCNTIAL FIT TC GAS TYPE2 



CC1220 

CC 1210 1=1,16 



0C1222 

121C CNT2{I) = £XPIC(lMC(2)*Rmi 



CC12 :3 

CALL LINPLT 1 R . CM2 , 1 6 , 1 , 0 , 1 2 . 1 , 0 ) 



CCl 2^ 2 

CALL CALPLT lC.,e..-3) 



C PLCT PATic OF GAS TVPEl TC GAS TYPC2 

CC12Af 

CALL ASCALE (R12,<i. , Id, 1, 10. ) 



CC1251 

CRG = SIG13)=R12I 17) 



CC125A 

5C=SIG14)=R12I1£) 



001256 

CALL AXES(0.,0.,0., 10., XCRG, X S , . 5 , 0. , EHCHANAF L S , . 1 5 ,- £ ) 



CC1271 

CALL AXES (0. ,C. ,90. ,A. ,CRG,SC ,.5 ,0. ,6FRAT ICS, .15,6 ) 



CC13C<i 

CALL LINPLT ( R , R 1 2 , 1 6 , 1 ,- 1 , 1 3 , 1 , 0 ) 



CC131A 

re 13 1=1,15 

• 


CC13U 

1 I'M 1 ) = T I^M2)=R m 



CC1321 

£IGll) = rtl2(I )-PRl2( 1 ) 



CC1323 

£IG12)=R121I ) + Efll2(I ) 



CC1325 

13 CALL LINPLT ( T I K , S I G , 2 , 1 , 0 , 0 , . 1 , 0 ) 



r FIT FCHIZCNTAL LINE TC RATIC CAT/ 

C01337 

5L‘' = WSUM = 0. 



CC13M 

re 1211 1=1,16 



0013A3 

1211 V>SU'' = WSUY + CR12{ I )**2 


- 

CCl 3A7 

L=0 


• 

C_C 1 3 tS 

re 1212 1 = 1,16 



C01352 

. . . i*L + l._ 




VO 



CC1254 


M2(L) = WSL^'/tR12m♦^^2 


CC1357 


va)-R12( I ) 

. 

CC12£C 


>(L)=R(I ) 


Uf;i26 2 

1212 

CCNTIKUE 

• 

CClid? 


L TCT=L 

• 

0C13tA 


V7S.''=R£S ID*SLN = C. 


CCl3t7 


CC 1213 L=1,LTCT 


CCl 3*fC 


HSX=VTSM + Vn12I L ) 


■' 'C013 /2 


5U»'=SU^'+Y(L)«V(12(L ) 


C C 1 2 ■< 5 

1212 

CCNTI NUE 


cri277 


Y I^TC = SUV/WTS^' 


v:ciAci 


tC 12U L = 1,LTC7 


CC1AC2 


9FSID=«ES in*W12(L)'M Y(L)-YIMC1**2 


C C 1 -V 1 2 

12U 

CCN7IMJE 


cuiAie 


ST=R = £0RT(RSSIC/( WTSf'fl LTGT-l ) ) ) 


C01425 


rC 1215 ! = 1,U 


CCH2« 

121E 

•<12(1) -YIMC 


C PLOT STRAIGHT LINE FIT 

CC1A21 


CALL LINPLT ( R t P 1 2 , 1 6 » 1 tO , C t 1 1 0 ) 


UCU^^ 1 

• 

CALL NOTATE ( 12 . 1 1 . C » .2 , 2 3HP A T I C STRAIGHT LIKE PIT»C.»23) 

. 

CCl*A5 


CALL KCTATE 112. 5, . 6 1 5, 1 2FY- 1 NT ER C EP T* , C 12 ) 


CCU5 1 


call NUVBER lU.l, .6,.15,YIKTr.O..A) 


001 AE 5 


call K'CTATE (12.5, .2, .15. 15HSTAN0A9C ERRrP= . C. . 1 51 


CCl^t 1 


call KUMRFR (U.5, .2,.15,STER,0. ,A) 


• coi^ies 

9*5 

IF (NCS2.EQ.0I CALL CALPLT (C.,8.,-2) 


ecu 71 


call NCTATE (12.,3.5,.?,AFCATE,0.,A) 


OOU75 


CALL NOTATE ( 1 2 . S , 3 . 5 r • 2 • DA T E , C . . 8 ) 


CCISCI 


CALL KCTATE ( 1 2 . , 3 . 1 , .2 , A FNC C E . Q . . A 1 


CC15CE 


CLAL^IUUAL 


CCl ;>C7 


CALL MJMOEP (12. 9, 3.1, .2, COAL.O..-1) 


GC1513 


IF (ISTL.EC.l) CALL KCTATE ( 1 2 . , 2 . 7 , . 2 . 2 1 J-S T G . LAMP CALIPPAT 

IQN.O. 

fill) 

CC1E21 


CALL CALPLT (C.,0.,999) 


CClf 2A 


STOP 


CC1E2E 



• 


VO 



CCCCIC 


CC0CA2 


CCCC46 


CCCC47 


OCCOf 1 


ccocet 


CCCC <1 


CC0073 


CCCCT4 


S L B R C UT I nTTSTn S ( F "C P S , F R F L t F ,'R E S C', N C S ) 


FCF(Z) = Z03S-Z’»EXP(Z/G) 


FCFP(Z)=EXP(Z/G)*(Z/F-l.) 
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